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Abstract 

We investigate the existence, uniqueness, and L^-contractivity of weak solutions to a porous medium 
equation with fractional diffusion on an evolving hypersurface. To settle the existence, we reformulate 
the equation as a local problem on a semi-infinite cylinder, regularise the porous medium nonlinearity 
and truncate the cylinder. Then we pass to the limit first in the truncation parameter and then in 
the nonlinearity, and the identification of limits is done using the theory of subdifferentials of convex 
functionals. 

In order to facilitate all of this, we begin by studying (in the setting of closed Riemannian manifolds and 
Sobolev spaces) the fractional Laplace-Beltrami operator which can be seen as the Dirichlet-to-Neumann 
map of a harmonic extension problem. A truncated harmonic extension problem will also be examined 
and convergence results to the solution of the harmonic extension will be given. For a technical reason, 
we will also consider some related extension problems on evolving hypersurfaces which will provide us 
with the minimal time regularity required on the harmonic extensions in order to properly formulate the 
moving domain problem. This functional analytic theory is of course independent of the fractional porous 
medium equation and will be of use generally in the analysis of fractional elliptic and parabolic problems 
on manifolds. 

Dedicated to Juan Luis Vazquez on the occasion of his 70*^ birthday 


1 Introduction 

For each t G [0,T], let T{t) C be a smooth and compact d-dimensional hypersurface without boundary 

evolving with a given velocity field w. In this paper, we are interested in the well-posedness of the fractional 
porous medium equation 

u{t) + (-Ar(t))^/q'u™(t)) -b u(t)Vr(t) • w{t) =0 on T{t) 

u(0) = uo on Fo := F(0) 

for m > 1, where uq is a given initial data, u™ := \u\^~^u as usual, and (—Ar(t))^^^ is the square root 
of the Laplace-Beltrami operator on F(t), which is a nonlocal first order elliptic pseudodifferential operator 
[48, 50, 61, 54]. 

If the fractional Laplacian in (1) is replaced with the ordinary Laplace-Beltrami operator — Ap((), (1) 
would be a porous medium equation on an evolving surface. Porous medium equations on stationary domains 
have, of course, attracted a considerable and well-developed literature. We refer the reader to the book [56] 
by Vazquez which is a comprehensive study of the mathematical analysis of the equation (and it also contains 
many references) and results on the porous medium equation on manifolds can be found in [56, §11.5] and 
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[12]. We will also say a few words about the non-fractional moving case in the conclusion of this paper. The 
investigation of fractional porous medium equations was instituted in [26] where the authors examined such 
a problem on involving the square root of the Laplacian and gave a complete theory of the equation, and 
indeed, our work is motivated by the results in that paper. There, the existence was proved by discretisation 
in time of a localised formulation of the equation and then the application of the Crandall-Liggett theorem 
[24]. Those results were generalised in [27] to a wider range of fractional powers of the Laplacian (—A)'* with 
exponent s G (0,1) on a stationary domain il C using the extension method introduced by Caffarelli and 
Silvestre in [19]. Existence was proved in [13] (for a more general nonlinearity) in a different way through the 
theory of semigroups and maximal monotone operators. Our model (1) differs from all of the aforementioned 
works since it is on a moving space. 

Other related works in the literature include variants of nonlocal porous medium equations such as those 
with variable density [47, 46] and different fractional operators [9]. We also mention [5, 20, 53, 43] where 
elliptic fractional problems are studied in the setting of the Laplacian on a bounded domain with Neumann 
boundary conditions, and [36] where a degenerate parabolic equation arising in crack dynamics is considered, 
again in the Neumann setting. One can also find numerical and finite element analysis for elliptic and parabolic 
problems in [44, 45]. As is evident, there has been an extraordinary amount of activity in fractional diffusion 
problems in the last decade or so. A good survey of recent and current output involving nonlinear fractional 
diffusion can be found in the articles [57, 58]. 

In terms of the analysis, a common preliminary step when working with half-Laplacians is to rewrite the 
problem locally using a Dirichlet-to-Neumann map [18, 7, 51, 23]. We will also reformulate (1) using such 
a map; this step is likewise performed in [26, 27] but from here on, the type of approaches taken in [26, 27] 
are problematic in our setting because of the additional complexity engendered by the evolving domain. 
For example, one could attempt to pull back the problem onto a reference domain (the resulting expression 
is not too cumbersome if the evolution of r(t) is prescribed particularly agreeably) and try to employ an 
appropriate time-dependent version of Crandall-Liggett [25, 31, 39] to the resulting equation (which will have 
time-dependent coefficients) but these theorems are difficult to apply even when the evolution of the domains 
is highly simplified. Therefore, we choose a different way to approach this problem, which we shall outline 
below, starting from the foundations. To our knowledge, the type of approach developed in this paper has 
not been used before in the fractional setting, even in the stationary case. The challenges and peculiarities 
that arise due to the moving domain will be highlighted in due course. 

Before we proceed, let us remark that fractional Laplace-Beltrami operators on various classes of manifolds 
have been studied in [7, 51, 23] through extension problems in the style of Caffarelli-Silvestre [19], but a 
convenient work detailing all the relevant properties of the half-Laplacian on closed manifolds in a Sobolev 
space setting appears lacking, so this paper is useful also in this respect. With this in mind, it is worth 
emphasising that the first part of this paper, comprising of §2-4, is independent of the second part which 
consists of §5 and § 6 , and indeed the reader can read the first part in isolation. The first part can be of use 
for other fractional diffusion problems on (evolving) manifolds and the second part can be thought of as an 
application of the first part. See the outline below for more details. 

1.1 Reformulation of the equation and main results 

A natural way to define (—Ar(t))^^^ is through a spectral definition which we describe now in greater generality. 
Indeed, suppose that 


(M, g) is a connected closed smooth Riemannian manifold 

and let ^k)ken be the normalised eigenpairs of the Laplacian —Am so that —AM^Pk = ^kPk 
fc; it follows that 0 = Aq < Ai < A 2 < ... ^ oo and ipo = [ 37 ^ Theorem 3.2.1]. The pk 

orthonormal basis of L'^{M) and are orthogonal in For smooth functions it, define 

00 

( Am) ^ u := A^ {u,(pk)L^[M)Pk- (2) 

k=l 


(Am) 

for each 
form an 
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The operator (—can be defined in a weaker sense through the action 

OO 

:= ^Xl/^{u,ipk)LHM){v,^k)L'^(M) 




which is sensible whenever u and v belong to the Hilbert space 


H{M) := <^ u G L\M) \ < oo 

I fc=i 


endowed with the inner product 


(3) 

(4) 


{u,v)h{M) '■= ( m , w ) l 2 ( m -) + (u, tpk) L'^(M){v, ipk) L'^(M)- 

k=l 

It is useful to have a Sobolev characterisation of the space H{M); in Lemma 2.7, we will see that 

HiM) = = ^2^(M) = (L2(M), 

i.e., H{M) is exactly the fractional Sobolev space (see [55, §7.2.2, §7.3.1, §7.4.5] for more details on 

the second and third equalities). In the later sections, we will be working on hypersurfaces so it is convenient 
for our purposes to introduce the Sobolev-Slobodeckii space (where T is a sufficiently smooth 

hypersurface) defined using the Gagliardo norm (see [4] and references therein): 

TTi/ 2 . 2 ('p) _ g I ;= |w(a:)p dcr(a:) + ^ dCT(a:)dcr(y) < ooj . 

Of course, this space is equivalent to iJ^/^(r) with an equivalence of norms (see [60, §1.4.2 and Theorem 5.2 of 
§1.5.1], [40, Theorem 7.7, Chapter 1], [40, Chapter 1, §15] and [34, §1.3.3]), but it is important to distinguish 
between these spaces when T = T(t) is time-dependent because the constants in the equivalence of norms will 
depend on t in an unknown way. 

The spectral definition of (—Am)^^^ in (2) is not particularly amenable to a convenient theory of weak 
solutions; however, there is a way to localise the fractional Laplacian (see [7, 51, 23]). With C := M x [0, oo) 
and g denoting the trivial product metric on C, consider the problem 

AgV = 0 on C, v\oc = u (5) 

where dC = M x {0}. Whenever u belongs to H{M), the equation has a unique weak solution v = £u, called 
the harmonic extension of u. This harmonic extension £u belongs in general not to H^{C) but to the larger 
space 

X{C) := where ||u||^(c) := \\Vgv\\%^c) + WTvWI^^m) for v G H\C) ( 6 ) 

with T: H^{C) H{M) denoting the trace map onto M x {0}, so that £: H{M) —>• X{C) (this type of space 
X(C) was first defined in a different setting by Stinga and Volzone in [53]). As we shall see in Lemma 2.6, 
the fractional Laplacian is recovered as a Dirichlet-to-Neumann map: 


{{ — AmY^'^U.v) h{MY ,H(M) — 


_ ,u) 

y-^ / H(M)*.H{M) 


where = (0, —1) is the outward normal to C. All of this will be laid out in detail in §2. 

Setting 'L(r) := and C{t) := r(<) x [0, oo), the above characterisation implies that one can rewrite 

(1) as 

u(t)-kM(t)Vr(t) • w(t)-h = 0 ondCit) 

v{t) =£t{'i'{u{t))) 
u{0) = uq on To 
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where Et is the map E with the manifold M chosen to be r{t) and i'{t) = (0, —1) is outward normal to C{t). 
Regarding the regularity of {r(t)}tg[o,T]i we will assume Assumption 3.1 on p. 14 and that 

there exists a constant Ai > 0 such that Ai(t) > Ai for all t G [0,T] (A^) 


where Afe(t) denotes the fc-th eigenvalue of — Ar(t); see Remark 3.2. A proper weak formulation of this 
problem requires the use of appropriate functional spaces. In [2, 3] (see also [4]), we defined generalisations 
of the Bochner spaces LP{0,T;Y) to handle functions defined on evolving spaces: given a family of Banach 
spaces Y = {I^(t)}tG[o,T]j a family of uniformly bounded linear homeomorphisms {(pt'. Yq —>■ l^(t)}tG[o,T] 
with uniformly bounded inverses Y{t) ^o}tG[o,T]i and t i—>■ measurable for all u G lb, we, 

generalising some work by Vierling [.59], defined the Banach spaces Ly as 


Ll.= 


{u : [0, T] ^ UtG[o.T] Y{t) X {t}, {u{t), t) 

{u G I eSSSUPig[0,T] ||M(t)||y(t) < 00} 


<()_(.)«(•) GLnO,T; lb)} 


endowed with the norm 




(/o^llwWIlyp)) ^ forpG[l,oo) 

^esssuptg[o,T] ||u(t)||y(t) for p = oo. 


for p G [1, oo) 
for p = oo 


(Note that we made an abuse of notation after the definition of the first space and identified u{t) = {u{t),t) 
with u{t).) The space W(y, Z) := {u G Ly \ u G A|} with ii the weak time or material derivative refers to an 
evolving space version of a Sobolev-Bochner space; this notion will be properly defined in §3 where we shall 
also make clear the assumptions on the evolution of the hypersurface T(t). This theory will allow us to define 
the following spaces (amongst others) after we make and check the relevant assumptions in §3. 

Space Ly Formed from {y(t)}tg[o^r] 

Li.(c) {L^cmtrnT] 

{^HC(0)}tG[0,T] 

^X{C) {^(^(0)}tG[O.T] 

In order to obtain measurability in time of f Et{'i’{u{t))) for u G L^^yj 2 (recall that each Et was defined 
individually at each moment in time as the harmonic extension on r(t)), we will consider in §4 the “Ax(c) 
harmonic extension” problem: given u G L^i/ 2,2 j find Eu = v € L'^x{C) ®och that 


Space Ly Formed from {y(t)}tg[o,T] 
Ll, {Lnr(t))}te[o,T] 


AgV = 0, Tu = u 


(7) 


holds with T: L‘^x{C) -^^ 1 / 2,2 the trace map. Then we will show that (EM)(f) = Etu{t) for almost all t, 

which gives the desired measurability. Of course, in the stationary setting, this issue of measurability would 
not arise and there would be no need to consider (7). Now we can think about what we mean by a weak 
solution. In what follows, given p G L^i/ 2 , 2 j we denote by Ep G an arbitrary extension of p that 

satisfies TL ?7 = p. 

Definition 1.1 (Weak solution). A weak solution of (P) is a function u G with E(4'(m)) G 

satisfying 

[ p{t)u{t) datdt + [ [ Vg,^t)Sti'i'{u{t)))Vg(^t){Ep){t) datdt = [ uop{0) dao 

r(t) Jo Jc{t) Jto 

for all p G W(IT^/^’^, L^) with p{T) = 0. Here, dat means the surface measure on T(t). 



From now on, for brevity, we will omit the measures in any integrals. We will prove the following theorem 
in §6, which is the main result of our paper. 

Theorem 1.2 (Well-posedness of the fractional porous medium equation). Under Assumption 3.1 and (A^), 
given uq G L°“(ro), there exists a unique weak solution u G L|[bo nL^_iy2,2 to (P) with E('I'(it)) G L'^x{c) 
the sense of Definition 1.1). Furthermore, we have the following properties: 
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1. Boundedness: for all t £ [0,T], u[t) £ L°“(r(t)). 

2. Conservation of mass: for all t £ [0,T], 




Mo- 


3. -contraction principle: if uqi and M 02 are two pairs of initial data in L°°{Tq), then the respective 
solutions Ml and M 2 satisfy 



{Ui{t) - M 2 (t)) + 


< 


/ (moI — Mo 2 )~' 
Jro 


for all t £ [0, T], 


An immediate consequence of the contraction principle is the following. 

Corollary 1.3 (L^-continuous dependence and comparison principle), //mqi and M 02 are two pairs of initial 
data in L°°{To), then the respective weak solutions ui and M 2 of Theorem 1.2 satisfy the -continuous 
dependence result 

[ \ui{t) - U 2 {t)\ < [ |moi-mo 2 | for allt €[Q,T]. 
dr(t) Jto 

If uoi < M 02 a.e., then ui{t) < U 2 (t) a.e. in r(t) for all t. 

Let us discuss how these results compare to those in the stationary case considered in [26, 27]. Theorem 
1.2 and its corollary correspond to parts i, ii, iv and v of Theorem 2.2 of [26] and to Theorem 7.2 of [27] in 
the half-Laplacian setting. In terms of the proof, our methods are quite different, as already discussed earlier. 
Let us sketch the proof now. 


1.2 Plan of the proof 

In order to solve (P) and prove Theorem 1.2, we will first approximate the nonlinearity 'L by well-behaved 
smooth approximations and seek to solve (P) with replaced by 'kfe. This directs us to study the 
non-degenerate problem 

upit) + up{t)Vr{t) ■ w(t) -b = 0 on dC(t) 

vp{t) =£tW{ui3{t))) 

M/3(0) = Mo on To 

where /3: ffi. —>■ R satisfies 

/3(0) = 0,/3 is C^(IR) (and Lipschitz) 

/3',(/3-')',(r')"ei°“(K), and (A^) 

there exist constants > 0 with 0 > Cp’ and 

To show well-posedness of (P/?) one could try a Galerkin method but a complication involving the unbounded 
cylinder C{t) arises due to the surface evolution, see Remark 5.4; this suggests truncating the cylinder C{t) in 
the unbounded direction. So we consider in §2.5 a truncated harmonic extension problem and show that its 
solution approximates the (untruncated) harmonic extension in some sense: given u £ H{M), with Sru = vn 
denoting the weak solution of 


AgM_R = 0 on Ch := M X [0,R], VR\Mx{o}='a, vr\mx{r}=0^ (8) 

we will show in §2.6 that 'S/gSnu —>■ 'Vg£u in L^{C) as R ^ 00 . As with £t, we define £R^t as £r with 
M = r(t) and CR{t) := r(t) x [0,i?], and consider the following problem as an approximation of (P^): 

M/3fi(^) + upR{t)Vr(t) ■ w(t) -b =0 on r(t) x {0} 

V/3R(t) = £R^t(P{u/3R(t))) ^ 

M/3ii(0) = Mo on Tq. 
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We can define the spaces ^'h^{Cr) truncated cylinder just like before, and consideration 

of an truncated harmonic extension” problem like (7) in §4 will lead to a map and show the 

measurability in time of Sa^t- We will use the Galerkin method to solve (P^r) in §5.1, see Remark 5.2 where 
we explain the choice of our Galerkin approximation; this requires emphasis due to a technical difficulty in 
the evolution-dependent projection operators associated to the Galerkin basis. Then we will pass to the limit 
in R in §5.2 in order to settle (P/ 3 ) and the following theorem will be proved. 

Theorem 1.4. Under Assumption 3.1, (^4^), and {Ap), given uq € L°°(ro), there exists a unique solution 
Up S to (P,a) with U/ 3 ( 0 ) = uq and¥.{f}{up)) € satisfying 





up{t)T]{t)\7r{t) • w(t) -h 


/ / V,p)ft(/3(u/3(t)))V,p)(i?r?)(t)=0 

Jo Jc{t) 


(9) 


for all r] € T^i/ 2 , 2 ; where the duality pairing is between W and Furthermore, mass 

is conserved and the -contraction principle holds for almost all t € [OjT]. 

With P chosen to be the regularisation 4/^, this theorem gives us a sequence {u/c}feGN where Uk G 
satisfies Ufc(O) = uq, E(4'fc(ufe)) G ^x(c)^ equation (9) with /3 replaced by 

'kfc and Up replaced by Uk- Then we pass to the limit in k using energy estimates and the identification of 
limits is handled with the theory of subdifferentials of convex functionals in §6 where the proof of Theorem 
1.2 is concluded. 

In [26, 27], the authors prove results for existence with integrable data too, as well as other properties 
besides, including regularity, smoothing effects and extinction of solutions. As the next step to our results, 
studying regularity in time would be natural (and useful) but it appears difficult in our setting. We comment 
on this in more detail in the conclusion. 


1.3 Outline 

It is clear that we need to properly study the harmonic extension maps St and Sr^, which we take care of 
in §2 in the general setting of closed Riemannian manifolds. In §3 we shall check that the spaces Ly listed 
above are well-defined and prove some preliminary functional analytic results. We then study the maps E 
and Efl in §4. After this, we tackle the non-degenerate problem (P/ 3 ) in §5 and then prove the main theorem 
in § 6 . We will finish with some concluding remarks in §7. Let us emphasise that §2 is useful more generally 
for fractional problems on closed manifolds and §3-5 are useful for fractional diffusion problems on (evolving) 
hypersurfaces. Only in §6 do we specialise to the porous medium equation. 

1.4 Notation 

We use the overline “ in different contexts. When applied to functions u, it means the spatial mean value: 
typically u = u or u = jrlpj /pp) u. When applied to symbols like E or f like E or S, the meaning 

usually is that the map with the overline is a linear extension, for example, f is a linear extension of f to a 
larger space. Symbols of the blackboard bold style like E refer to maps between the evolving Bochner spaces 
Ly, whilst symbols of the calliographic style like S refer to maps between Sobolev spaces of the form 
The notation j-j denotes a seminorm; usually the part of the corresponding norm is omitted. 

As a convenience for the reader, we give here a list of the major notations and symbols that we use in this 
paper along with the page number of definition or first usage. 


Notation 

Page 

Notation 

Page 

(Am) 

p. 2 

St, Srp, St, Srp 

p. 16 

H{M) 

p. 3 

'Tit 'Tt 

p. 16 

X{C) 

p. 3 

T, T 

p. 16 

(A,) 

p. 4 


p. 16 

S, £ 

p. 8 

E, E 

p. 19 

Sr, Sr 

p. 11 

E/{, Efl 

p. 21 

^R 

p. 13 


p .21 
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2 The fractional Laplacian on compact Riemannian manifolds 

Throughout this section, we assume that {M,g) is a Riemannian manifold as given in [Am)- One aim of this 
section is to realise the fractional Laplacian on a closed Riemannian manifold as the Dirichlet-to-Neumann map 
of a harmonic extension problem in a Sobolev space setting. We will define an operator 8 : H[M) X[C) for 
this purpose. We also study the truncated harmonic extension by means of an operator Sr : H(M) — (Cr), 
and then prove that Eji approximates £. 

Remark 2.1. We do not consider the case where M is an open manifold (i.e., a manifold with boundary). 
If dM 0 and we place Neumann boundary conditions then most of what follows should be similar, since the 
eigenvalues and eigenfunctions behave similarly to the closed setting. If instead Dirichlet boundary conditions 
are taken then an analogue of the following results will hold; in particular one probably would not need to 
worry about differentiating between functions with mean value zero and those without, and the space H[M) = 
(L^(M), Wq’^(M))i /2 will be rather different. 

We will often be integrating or manipulating infinite series of functions term by term which can be justified 
by Abel’s test or the Weierstrass M-test. More details of this and lengthier calculations of what follows can 
be found in [1]. First, we begin with a brief discussion of Sobolev spaces on (semi-infinite) cylinders. 

2.1 Sobolev spaces on semi-infinite cylinders 

We can use the space II^[C) (utilised already in the introduction) defined in [6] as the linear subspace of Ll.iC) 
consisting of all v such that v and belong to lf[C), and endowed with the natural norm. Equivalently, it 
can be defined as the linear subspace of L^(0, oo; H^[M)) consisting of all v such that Vy S L^(0, oo; L'^[M)). 
This is precisely the type of Sobolev-Bochner space whose theory was developed by Lions and Magenes 
[40, Chapter 1, §2.2]. There is a bounded linear surjective trace operator T: II^[C) —^ H^^‘^[M x {0}) [6, 
Theorem 18.1], [40, Theorem 3.2, Chapter 1], possessing a continuous right inverse. Similarly, the spaces 
H^[Cr) can be dehned on the truncated cylinder Cr = M x [0,i?]. Theorem 3.1 of [40, Chapter 1] gives 
II\Cr) ^ (^^([O, i?]; iL^/^(M)), so that the linear trace operators TR,y=o, TR,y=R: II^{Cr) —?> il^/^(M) 
defined by [TR,y=ov)[-) := u(-,0) and {TR.y=Rv){-) := v[-,R) are also bounded. Furthermore TR^y=o is 
surjective [40, Theorem 3.2, Chapter 1]. 

Lemma 2.2. Ifv S II^[C), then y >->■ v{y) = v{y) is an element of II^[0, oo) and thus v S C°([0, oo)). 

Proof. A calculation verifies that v G ff^(0, oo), and Theorem 8.2 in [17] proves that each function in iL^(0, oo) 
has a unique continuous representative in C'°([0, oo)). □ 

2.2 Fractional Sobolev spaces and the fractional Laplacian 

The setting of a closed manifold is similar to the setting of Neumann boundary conditions on a bounded 
domain (see [53, 43, 52]), and now we motivate the definition of the half-Laplacian like [53, §2]. As mentioned 
in the introduction, let (Afc, (pk) be the normalised eigenelements of the Laplace-Beltrami operator —Am- For 
k 0, since [pk,Po)L^(M) = 0, ^ = 0. We also have = 1 + which implies that 

{ OO 

u G L^{M) 1 = ^(1 -f Xk)\{u,Pk)L^(M)\^ < oo 

k=0 


and for u G II^[M), one has 


OO 

-Amu = ^ Xk{u,(pk)L^(M)^k in 


with 

OO 

{ — ^MU,v)}j-^(M).m(M) = ''^Xk{u,ipk)L'^(M){'V,Pk)L'^{M)- ( 10 ) 

fe=l 
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With the Hilbert space H{M) as in (4), the previous two identities inspire us to define (—: H{M) —>■ 
H{M)* by (2) with the action (3). For u, v G H{M), it is easy to see the integration by parts formula 


Jm 

where ((-Am)^/'‘u, u) := ^l^^iu^^k)L^{M)iv,V^k)L^{M), and we have 

2 


fe=l 


2.3 The harmonic extension problem 

Recall the problem (5). If u = 1, then its harmonic extension is u = 1, so u i—^ u does not map into H^{C) in 
general. Therefore, we will work in the bigger space X(C), defined in (6). 

Remark 2.3. The constant functions belong to X{C). To see this, take c G R and Cn G H^{C) with 

(c :y G {0,n] 

Cnix,y) = lf^(2n-y) :yG{n,2n] 

[o :y G (2n,oo) 

which satisfies X mCu = 0 and dyCn = —c/nx(n, 2 n)iy)- Note that 



so c := (cn) is a Cauchy sequence in the X{C) norm, and it follows that 


ll^^llxcc) := llcnll^(c) = + c^ = Jim \M\ = \M\c^. 

Then c can be identified with the constant c. 

Lemma 2.4 (Extension of the gradient to X{C)). The gradient V^: H^{C) —?> L^{C) extends to a bounded 
linear map Xg : X{C) —?> L'^{C) which satisfies Xg\ffi(^c) = Ng and XgV = lim„_>oo XgVn for G H^{C) such 
that Vn ^ V in X{C). 

Proof. Clearly Xg: H^{C) —^ L^(C) satisfies ||Vgu||^ 2 (c) — Ihllx(c) ^ ^ N^{C). Since H^{C) is dense in 

X{C), the bounded linear transformation (BLT) theorem provides the result. □ 

Theorem 2.5. For every u G H{M), there exists a unique weak solution £u = v G X(C) to the harmonic 
extension problem (5) satisfying (£u)(-,0) = u(-) in Lf{M) and 


Xg{£u)XgTj = 0 for all rj G H^{C) with Trj = 0. 


When u = 0, we write the solution as £u which is such that 0 y ^ [0,oo). The 

map £: H(M) —^ X{C) satisfies £u = £(u — u)+u and Xg{£u) = Xg£{u — u). Furthermore (if u = 0), 
£u G C'°([0, oo); L‘^{M)) n (^“((0, oofH^M)) and 


IIP ||2 ^ II^IIl2(m) 


2A 


1/2 




L^{M) 


II^S^^II L2(C) 

Finally, the harmonic extension £u (for u = 0) is the unique minimiser of the energy 

J:{vG H\C) \Tv = u}^R, Jiv) := ^ J \Xgv\^ 


( 11 ) 

(12) 
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Proof. The proof of the well-posedness is essentially the same as that of Theorem 2.1 in [53]. Suppose for 
now that u = 0. Set 

OO 

» ^ \ 1 / 2 

{£u){y) := v{y) := ^ e~^ (u, <fk)L^(M)Pk, 

k^l 

which is a sum that converges in L^(M) for each fixed y € [0,oo). It satisfies 

/ iv(y)j^ l(u,(pk)L^(M)j^ and f \V gv{y)\'^ = 2^ \{u, ipk)L^(M)\‘^ 

^ k>l k>l 

where we used (10). These expressions can be integrated over y term by term since the sums converge 
uniformly, and doing so leads to properties ( 11 ) and ( 12 ). 

Now let rj G H^{C) with Trj = 0. For almost all y, r]{y) = Pk)L^{M)Pk] let us write the 

coefficients as r]k{y)- We see that 


f °° 1/2 

/ ^Mv{y)VMviy) = ^Xke~y^>‘ iu,ipk)L^M)Vk{y) 

k=i 


= / 

Jm 


iiy)viy) 


using (10). By integrating by parts we get 


/ Vmv'VmV + VyTjy = / y mvS/mV -'k^vvV + / VyV = 0 
Jc Jc JdC 

as V has zero trace. This proves that ?; is a weak solution. Uniqueness follows by taking the difference of the 
weak formulations satisfied by two solutions and testing with the difference of the two solutions (which has 
trace zero). Therefore, the map £: {u G H{M) | m = 0} —>■ H^{C) is well-defined. Now suppose u 0. Define 

£{u) := £{u — u) + u. 

Note that Vg{£u) = Vg£{u — u)+ Vgll by linearity and the fact that £{u — u) G H^{C). Let us choose 
Un = Cn G H^{C) as in Remark 2.3, which tells us that lim„_>oo ^gUn = linin-s-oo ~^X(n, 2 n)(y) = 0 in L'^{C), 
i.e., VgU = 0. This proves that Vg{£u) = Vg£{u — u). 

For the minimisation property, take w G H^{C) with Tw = u, test the weak form £u = v satisfies with 
V = V — w and use Young’s inequality: 

II^s^IIl2(c) — 2 + 2 


and rearranging shows f£(v) < J{w). Uniqueness follows since J is strictly convex. 


□ 


We will often (but not always) write Vg instead of Vg. From (12), we find that £: H{M) — X{C) is an 
isometry: 

The next lemma is fundamental (see also [18, 51]). 

Lemma 2.6. The fractional Laplacian of u G H{M) is recovered through the Dirichlet-to-Neumann map: 

{—Am)^^^u=— lim in H{M)*. 

y_>o+ dy 

Proof. If u = 0 and v € H{M), taking the limit ?/ —>■ 0+ (using Abel’s test) in 

OO 

E l /2 

V '' {u,Tk)L^{M){V,Tk)L^{M) 

fe=l 

and comparing the result to (3) gives us what we expected. The case u 7 ^ 0 follows easily. □ 
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Lemma 2.7. The space H{M) = with an equivalence of norms. 

Proof. Given u G H{M) with u = 0, define v = 8u, which we know belongs to H^{C) from Theorem 2.5 and so 
Tv = u G since T has range in For the case it 0, we have that u = u — u+u G 

Now we prove the reverse inclusion. Recall that a function u G Lf{M) + H^{M) belongs to the interpolation 
space as defined by the if-method if the following norm is finite: 




(M) - 




1/2 


uoeL^(M) 

uieH^iM) 


See [40, Chapter 1, §15], [16, Appendix B], [41, Appendix B] for more information. We follow the ideas of the 
proof of Theorem B.2 in [16] now. Let u = '^kPk G and v = G H^{M). Then 




^ Mhhm)) 


= inf ^ itfc - itfc +t ^(1 + Afc) itfc 
^ ' \k=0 fe =0 


t^(l+A/c) 2 

= 2^ ^|ufc| 


fe=0 


1 + t^(l + Afe) 


because the expression in the infimum is minimised when Vk = itfe/(l + t^(l + A^)). Therefore, 


LXJ ^OO ^ ^ 

ll“llili/2(M) = ■*" ^k)\uk\'^ X + 1^(1 + Afc)*^^ ~ '2 ^ \/l + Afejltfe]^ > — lltl//(M) 

k —0 /c—0 

which implies that l|itl|//(M) ^ 7 r “^/^(2 + 7 r)^/^ \W\\m/^(M)- above calculation, using + Afe < l + ^/\k 

shows that 1 |m 11 jji/ 2 (m) < 2 - 1 / 27 ^ 1/2 . □ 


We could also have proved this lemma via the J-method of interpolation [41, Appendix B] and Weyl’s 
law [-37, Chapter 3, equation (3.2.24)], as is done in [13, §3.1.3] on a bounded domain. Another approach, 
relying explicitly on the Gagliardo norm on i7i/2(M) when M is a hypersurface, can also work with the use of 
two-sided Gaussian estimates on the heat kernel, similar to [53, §2.2] for the case of the Neumann Laplacian 
on a bounded domain. 

We introduce the following cut-off function which will be useful here and in § 6 . 

Definition 2.8 (Cut-off function). For any p > 0, there exists a smooth cut-off function ifp such that 


^p{y) 


1 :yG[0,p] 

0 :yG [ 2 p,oo) 


and —^Y^'0(1 — j) < tp'piy) < 0 on [p, 2p], with C not depending on p. 

Define a map M: iJ-i/2(M) by 

(A/m,/ i)^-i/2(7Vf),_f/i/2(M) := y Tg{£u)\7gh (13) 

where h G H^{C) is any extension of h (i.e., Th = h). This map is well-defined since if we had two arbitrary 
extensions hi and / 12 , then 


J Vg{£u)Vghi- j Vg{£u)Vgh2= J Vg{£u)Vg{hi - h 2 ) = 0 

by definition of £u and since T{hi — /i 2 ) = 0. The fact that the extension can be arbitrary will be extremely 
useful later on. Furthermore, by choosing in (13) h = £{h — h) -\- fjph G H^{C), one can see that A/m is 
linear and that it is indeed in the dual space of i7^/^(M). We can write A/m = d£u/diy\^^^, i.e.. A/ is the 
Dirichlet-to-Neumann map; this notation is justified since, if for example Ag£u G Lf{C), the standard Green’s 
formula implies wd£uldv = J^wAg£u-\- j^Vg£uVw = j^Vg£uVw. 
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2.4 Trace maps 

The trace map can be extended to the space X(C). 

Lemma 2.9. There exists a bounded linear trace map T : X{C) —>■ H{M) such that 

< ||r;||x(c) for v e X{C), 

Tw = Tw if w € H^{C) C X[C), and Tw := lim„_>ooT?«n for Wn S H^{C) converging to w in X{C). 

Proof. This is similar to [53, Lemma 2.4]. Let w G H^{C) be arbitrary with Tw =: wq. If rco = 0, by (12), 

kolff(M) = < 2J{w) = ||Vga;||^2(c) , 

and since this inequality involves seminorms, we can drop the assumption wq = 0. Adding ||wo||^ 2 (^) to 
both sides shows that T: H^{C) —?> H[M) satisfies \\Tw\\uf^f^~^ < ||rc||^(c) ■ Then the BLT theorem gives the 
result. □ 

The following lemma is a seminorm boundedness property of the trace map; note the Gagliardo seminorm 
on the left hand side (the proof of Lemma 2.9 had the H{M) seminorm on the left hand side instead). 

Lemma 2.10. Let M = T he a hypersurface of class . For every v G H^(C), 

|Tu|^^l/2,2(r) < C II VgUjl ■ 

Proof. If w € H^{C) satisfies v{y) = 0 for all y, then using the trace theorem and Poincare’s inequality on T, 

l|T'u||jyi/2,2(p) < Cl ^11 Vrt'||p2(p) + < C 2 l|Vg^lli2(p) ■ 

Now suppose that v G H^{C) does have not have spatial mean value zero for a.a. y. Then define v{x,y) = 
v{x, y) — v{y) which satisfies Jp v{y) = 0 and v G H^{C) by Lemma 2.2. Then, by the above inequality. 


l|T^||^p^l/2,2(•p) < C 2 Ml + |p| J dyv{y) 


L^{C) 


but, using Lemma 2.2, the left hand side is greater than |Tz^|ppi/ 2 , 2 (p) 
constants. 


^ <G3||V,u|M2(c), 

because the seminorm does not see 

□ 


2.5 The truncated harmonic extension problem 

Define Hq{Cr) := {p G iJ^(Ci{) | TR,y=o?7 = TR,y=RP = 0}; this is a Hilbert space because it is a closed linear 
subspace of H^{Cr). For p G Hq{Cr), it follows by Fubini-Tonelli that for almost all x, p{x,-) G IIq( 0,P). 
Thus the Poincare inequality on [0, i?] implies for almost all x that 

rM(x)p<c'p ripy(x)i^. 

Jo Jo 

Using this fact in the definition of the norm of Hq{Cr) gives IMIIpi(Cji) — ^ so that l|Vg-|M 2 (p^) 

is an equivalent norm on Hq{Cr). 

Theorem 2.11. For every u G H{M), there exists a unique weak solution £ru = v G H^(Cr) to the truncated 
harmonic extension problem (8) satisfying (£ru)(-, 0) = u(-) and {£ru){-, R) = 0 in Lf{M) and 


' Cr 


VgVVgP = 0 


for all p G Hq{Cr). 
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When u = 0, we write the solution as Eru which is such that Jj^{£Ru){y) = 0 for all y € [0, i?]. The 

map Er: H{M) —>■ H^{Cr) satisfies Eru = £r{u — u) + ^^u. Furthermore, Eru G C°([0,i?];L^(M)) fl 
C°((0,E];B^(M)), dyERuG C°((0,B];L^(M)) and 

- ~~Tj2 II«-“IIl2(m) +4-R|M||u|2, (14) 

l|V§^fi«irL^(CH) ^ (l + l/(2sinh2(A}/"i?))) + (15) 

Proof. Suppose that u = 0 and define 

OO 

{ERu){y) := v{y) := ^ {ai{k,R)e'^y + cx 2 {k,R)e~'^y^ {u,Tk)L^{M)Tk 


where 


ai{k, R) = -- 


3 y/Tj^R 


and a 2 (k, R) = 


spTk R 


Qy/FkR — f^ — yJTZR ’ (^y/FkR — ^—y/\kR 

The formula for Eru comes from separation of variables and the infinite sum converges in L‘^{M) for all 
y G [0, i?]. We have 


/ k(y)P = / '^{ai{k,R)e'^y + a2{k,R)e~'^A \{u,ipk)L'^{M)\‘ 
JCn Jo ^ J 

~ o rr- „2,/ArR , „-2,/xrR o Iw^TkjL^^M)] 


k= 


2\/\f gZVFkR -(- g 2y>TR _ 2 


< 


2^/Xi 




and 


IM 


OO 

ju(y)p = 2^ Afc (ai{k, Rfe'^'^y + a2ik, \{u,y:>k)L^(M)\‘^, 


fc=l 


and formally, integration over [0, i?] of the latter quantity yields 

^2y/X^R _ —2y/X^R 

ivj.f = >; x/Ai-' 

^ Cr 


r ^ ,_ p 2 V^kR _ p-^V^kR 

Jcn ^ + e-2VWiJ _ 3 ^ 

“ _ g2v^-R _|_ g-2v^_R 

< X, p2V>TR + e-2V^R _ 2 ^ I 

OO 

X ‘/^fe)L2(M)P- 


fc=l 

< (1 


2sinh^(vVi?) 


k=l 


(16) 


(17) 


To make this rigorous, we should have integrated over [e,i?] instead of [0,i?] and then passed to the limit as 
e —^ 0 using the monotone convergence theorem on the left hand side. 

To see that is a weak solution, take a test function rj G Hq{Cr) with -qljj) = YAk=oAl{y)^Tk)Tk and 
calculate (using (10)) 


/M 


Vr'fVr?? = X] 

fc=i 


\k{ai{k,R)e'^y + a2{k,R)e {v{y),(pk)mM)iu,‘Pk)mM) = [ VyyV 

^ ^ Jm 


Then 



Vrr'Vr?? 



Vyyp = - 



VyPy + 



VyPy 
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with the last equality since 77 vanishes on the boundary; this implies the result. For u with rt 7 ^ 0, we set 
£ru := £r{u — u) + This is a solution because 


IM 


R-y -^ 

R “ 





u d 

R Jo dy 


f y = -ll {v{R)-vm = o- 

Jm ^ JM 


If u 7 ^ 0, (14) and (15) follow from (16) and (17) by noting that £r{u — u) _L (i? — y)u/R in L'^{M) and 
Vg£R(u — u) _L Vg(i? — y)u/R in L‘^{M) respectively, pointwise in y. □ 

Remark 2.12. Define a form qr: H{M) x H{M) —^ R 6?/ 


aR{u,f]) 



VgSRUVgfj 


where fj S Hq{Cr) is an (arbitrary) extension of rj; the choice of extension does not matter, since for any two 
such extensions 7)1 and 7 ) 2 , 



VgSRUVgfjl 



Vg£RUVgf]2 



Vg£RUVg{f]l 


m) = 0 


by definition of the weak solution, because 7)1 — 7)2 S Hq{Cr). 


2.6 Decay and convergence of solutions of the truncated problem 

In order to compare functions defined on Cr and C, we define the zero extension 


^R-.{y&H\Q,R)\y{R)=Q}^H\Q,<x) by {2£Ry){y) = ^ (18) 

I 0 : otherwise 

which is an isometry. Clearly, we can also view ^r as a map 3£r : {77 G H^{Cr) \ 77 ( 37 , i?) = 0} —^ tJ^(C) and 
this is also an isometry. 

Lemma 2.13. For all u € H{M), 

\\Vg(£u - flu)||^ 2 (c) < \u - u\h^m) + 

Hence ^r£ru —>■ £u in X(C) as R ^ 00 . 

Proof. Firstly, let tjr = {£u — £ru) — £u{R)y/R which satisfies f]R{0) = r]R{R) = 0, and consider the difference 
of the weak formulations of £ru tested with rjR and £u tested with t^RrjR: 

0 = / \yg{£u-£RU)\'^-\7g{£u-£RU)Vg 

Jcr 


^u{R)y \^ 


so that 


'Cr 


\Xg{£u - £ru)\^ < / \Xg{£u - £Ru)\\Vg{£u{R))\ + \dy{£u - £ru)\\£u{R)\ 


'Cr 

< J \\Xg{£u - £_ru)"^ 


R 


L^iCn) 


'Cr 


\Vr£uiR)\^ +—\£u{R)\ 


where we used ab < a^/4 + 6^ on both products. Now, recalling that £u{R) = X)fe>i ® ^'^{u — u,ipk)ipk + u, 


Hr 


\£u[R)\^ < Re 


— 2i?\/Ai 


II«-“IIl=(m) + R\m\\u\' 
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and 


f \Vr£u{R)\'^ = ^ RXkC - u,ipk)\'^ <^y^e - u,ipk)\'^ < e 

k>l k>l 

hold (where we used using xe~^^ < e““), giving 



£Ru)\^<2e 


^ -2flVAr 

R 


11^ ^IIl2(M) 


2|M||?I|2 

R 


Secondly, note that 




poo _ 

= 2 Afce~^^^|(M - u, (/?fc)p 

fc>i 


_ g—S-R- n/XT 


1^ ^Ir(M) ■ 


Adding this and (19) implies the result. 
Lemma 2.14. For all u G H{M) with u = 0, 


||Sk£’RM-fM||i2(c) < Cp ^3e ^'^|m|p(m) + ;|e \\u''^ 


L^{M) 


^—2Ry/Xl ^ 


‘2'/Xi 


^IIl2(M) 


(where Cp is the Poincare constant on M). Hence SFrEru —>■ £u in L^{C). 


Proof. If u = 0, then £u{y) = £Ru{y) = 0 for all y. Therefore, with Poincare’s inequality on M, 


pit p pit p poo p 

/ / \£u-£ru\^<Cp / \Vg£u-Vg£Ru\'^<Cp / \Vg£u-Vg£FR£Ru\'^. 

Jo Jm Jo Jm Jo Jm 

Over the interval (i?,oo), we have 


( 19 ) 


□ 


r°° r r°° r p- 2 R\/Ar -2rF7F 

/ / \2fR£RU - £u\^ = / \£u\^ = Y^ \{u,(fik)L^iM)\^ < ^ \\u\\l^m)- 

Jr Jm Jr Jm 

Adding these two estimates and using the previous lemma yields the result. □ 

The next lemma describes continuous convergence. 

Lemma 2.15. If ur, u G H{M) with ur ^ u in L‘^{M) with ur = u = 0, then J!Fr£rUr —>■ £u in L^{C). 
Proof. Writing 3£r£rur — £u = 3fR£R{uR — u) + 2fR£RU — £u, using the triangle inequality and (14), 

ll-^R^'R'aR ~ ^r|| L'^{C) — ^ ~ ^11 L^{M) \\^r£rU — £u\\ p2('c) 

which tends to zero by Lemma 2.14. □ 


3 Function spaces on evolving hypersurfaces and preliminary re¬ 
sults 

We start with conditions on the prescribed evolution, in addition to {A\). 

Assumption 3.1. For each t G [0,T], let r(t) C be a smooth and compact d-dimensional hypersurface 

without boundary, and assume the existence of a flow $: [0, T] x —>• R'^+^ such that for all t G [0, T], with 

To := r(0), the map $j(-) := $(t, •): Tq —!> r(t) is a -diffeomorphism that satisfies (')) 

and 4?q(-) = Id{-) for a given velocity field w: [0,T] x R'^^^ —>• R'^+^j which we assume satisfies the 
uniform bound |Vr(t) • w(t)| < C for all t G [OiT]. A normal vector field on the hypersurfaces is denoted 
by [0,r] X R^+i ^ R'^+h 
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It follows that the Jacobian := det is and uniformly bounded away from zero and infinity. We 
denote by <i>g: r(t) —7> Fg the inverse of and define |r| := max^gfo^r] |r(t)|. 

Remark 3.2. The assumption (^a) is satisfied if for example each r(t) has non-negative Ricci curvature, 
or if the Ricci curvature ofT{t) is greater than p(t) < 0, where —p(t) < p holds for all t G [0,T] with p a 
constant. See Theorem f.d.l in [37] and the discussion afterwards. Also, instead of assuming (Aa), one could 
study the possible continuity o/1 i—^ Ai {t) through the theory of perturbations of linear operators [38]. Let us 
furthermore remark that all functional analytic results in this section not involving the harmonic extension 
maps remain true for r(t) of class . 


3.1 Function spaces 

In order to define the spaces Ly mentioned in the introduction, we need simply to verify a few assumptions. 

3.1.1 Spaces on the surface F 

For u: Fg —7> M, define {(j)tu){x) := {(j)r,tu){x) := u($q(x)). Fortunately, we already checked that the spaces 
and F^i/ 2,2 are well-defined in [2, §2.2.1] and §4.1 and §5.4] respectively. Recall from [3] that 
u G F^i/ 2,2 is said to have a weak material derivative ii G F^_i/ 2,2 if 


fT i-T 

C(i))w-i/2,2(r(t)),vvi/2,2(r(t)) = - 


f / CitMt)- / wWCWVrp) • w(t) 
0 Jr{t) Jo Jr(t) 


holds for all (] G {C S F^i/ 2,2 I '/>r,-(.)C(‘) £ ^^((0:1’); fF^/^’^(Fo))}, where ([ belonging to this set has a 
strong material derivative defined by C,{t) := (j)r,t{d/dt{(j)r-tC{t))) (see also [4]). In [4, §5.4.1] the evolving 
Sobolev-Bochner space 

W(VFi/2.2,ptc-i/2.2) - {y g ^ I y g L^_i/2,2} 

was shown to be well-defined and isomorphic (via 0r,-(.)) with an equivalence of norms to 
>V(1F1/2.2^ tR-1/2.2) _ I - g 2,2(0^ T; \ h' G ^^(0, T; 

and this implies that W(IF^/^’^, A L'^ 2 . 

The following lemma (which is surprisingly non-trivial) is useful later on; the proof of the continuity is 
the same as in Lemma 2.5 of [14] with the obvious modifications. 

Lemma 3.3. For a sufficiently smooth hypersurface F, */ /3: K. —>■ R is Lipschitz with /3(0) = 0, then the map 
fi: 1F^/^’^(F) —>■ 1F^/^’^(F) defined by composition is (sequentially) continuous and satisfies 

ll/3(M)llwi/2,2(r) < Lip(/3) l|M|lwi/ 2 , 2 (r) for all u G VF^/2’^(F). 

3.1.2 Spaces on the cylinders C and Cr 

Recall from the introduction that C{f) = F(f) x [0,oo), and set Cg := C(0). Given v G L^(Cg), define 
{(j)c,tv){x,y) := u(4>(,(x),?/). We have 




LHC{t)) 


/r(i) 


K4>‘(x),2/)P = 




To 


( 20 ) 


so fic.t - Lfi{Co) Lfi{C{t)). The inverse mapping is fic.-t - Lfi{C{t)) L^{Co) given by {4>c,-tw){x,y) = 
w{^'[{x),y) and these maps are linear homeomorphisms. Also, we see that if n G H^{Co), 


\fic,tv\ 


m(c(t)) 


lr(t) 


\37-,vmx),y)\^ 


lr(t) 


{x){Vroviy) o (x)) -f dyv{<i>l{x), y)]" 


<C| 


To 


TTCo) 


1(D4>(,)T o 4>0(z)(Vrof ( 2 :, y)) + dyv{z, y)\^J^ 
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( 21 ) 

( 22 ) 


Overall, we have shown that H^{Co) —>■ H^{C{t)) is bounded uniformly and well-defined. Finally, we 
have from (20) and (21) that t !->■ \\4’C,tv\\%i(^c(t)) continuous. Using the theory in [3, §2.2], this allows us to 
define the spaces U| 2 (c) (just ignore the gradient term). Clearly, the same argument allows us to 

define using a map 4>CR,t defined in the same way. 

Definition 3.4. We denote by £t and the maps £ and £ii defined in Theorems 2.5 and 2.11 respectively 
with the manifold M chosen to be Tft) (and likewise without the overlines). Similarly, we denote by Tt, 
TR,t,y=o and TR,t,y=R the trace maps T, TR,y=o and TR,y=R defined in Lemma 2.9 and in ^2.1 respectively 
with the choice M = r(f). 

Lemma 3.5 (Commutativity of the trace and pushforward maps). The following identity holds: 

Tt{(j)c,tv) = frA'^ov) for all v G H^{Co). 

Proof. We have ((c tV = v o £ H^{C{t)) and so TAc tV = u($q(-), 0), whilst on the other hand, tTov = 
u(-,0)oci>‘(.) = u(4>‘(-),0). ’ ’ □ 

Lemma 3.5 implies that if u G H^{C{t)), then, using the boundedness of 

||7tu||^yl/2,2(r(t)) < Cl \\'To4>C,-tV\\y^l/2,2(^Yo) — ^2 ll<('C,-t'c||^l(CQ) ^ C 3 

because of the trace theorem and the equivalence of norms between i?^/^(ro) and lU^/^’^(ro). This shows that 
Tt: H^{C{t)) —7> iU^/^’^(r(f)) is bounded independently of t. By the same argument, the uniform boundedness 
of TR,t,y= 0 : TR,t,y=R '■ H^{CR{t)) —>■ (r(t) ) also holds, and a version of Lemma 3.5 holds for these maps, 

which allows us to define Now, by Lemma 3.5 and (22), for v G H^{Co), we have 

\\^c,tv\\x(^cit)) ~ + l|7t^c,tu||j^2(r(t)) < C" 

which shows that (j>c,t' H^{Co) T{C{t)) has a useful boundedness property which, by the BLT theorem, 
allows us to extend fc.t to a bounded linear map (l)c,t ■ X{Co) —^ X{C{t)) defined as 

^C,tXo := lim fc.tVn in X{C{t)) for u„ G H^(Cq) with xq in X{Co). 

n—>-oo 

We also have the measurability of f 1 —>■ = linin->.oo \\4’C,tVn\\x{c(t))- Thus is also well- 

defined. Similar arguments can be made for the inverse operator of 4>c,t, denoted 4>c,-t' X(Clf)) —^ X(Cq). 
By a density argument, exploiting the continuity of the operators involved, we can show that 

Ttifc.tv) = fr.tiTov) for all u G X(Co). (23) 

3.1.3 Superposition trace maps 

Lemma 3.6. There exist bounded linear trace operators T: L^i/ 2,2 and T: Lx(c) T^i/ 2,2 

satisfying (Tu)(t) = Ttvff) and (Tu)(t) = Ttvlf) for almost every t. 

Proof. For v G L'^xiC)' tioiiue (Tu)(f) = Ttv{t). Then (Tu)(t) = Ttvff) = 4>r,tTo^c,-tv{t) by (23) which gives 
measurability in time, and we have the bound 

||(T^^)(^)|lM/i/2,2(r(t)) < Cl ||ro^c.-tu(t)||^i/2,2(r„) < C 3 \\fc,-tvit)\\^^^^^ < Q ||u(<)||^(c(t)) 

where for the second inequality we used the equivalence of norms between VF^/^’^(ro) and H{To) and Lemma 
2.9. This proves that T: ^^ 1 ^ 2,2 is well-defined as a bounded linear operator. The operator T can 

be seen as the restriction of T to □ 

By the same reasoning as above, we can prove the following lemma. 

Lemma 3.7. There exist bounded linear trace operators TR^y=o,TR^y^R: —>■ L^i/ 2,2 satisfying 

O^R,y=ov)it) = TR,t,y=ov{t) and {TR,y=Rv){t) = TR,t,y=Rv{t) for almost every t. 
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3.1.4 Some uniform bounds 

When we work with a time-dependent manifold M = r(t), we would like the constants in the gradient 
bounds (12) and (15) to be independent of time. The space i/^/^(r(t)) is equivalent to lT^/^’^(r(t)) with 
an equivalence of norms, as we mentioned in the introduction. However, the constants in the equivalence of 
norms result will depend on t and we have no information as to in what way the dependence is. This means 
that one has to be careful whenever one uses estimates from §2 involving the or norm in 

the evolving set-up. For this reason, we need the bounds in the next two lemmas. 

Lemma 3.8. There exists a constant C > 0 such that for all t G [0,T] and all u G H^^‘^{T{t)), 

II Vg(t)£^tu||i 2 (c(t)) < C* ll« - ■“llwi/ 2 , 2 (r(t)) ■ 

Proof. Let u = 0 and set U{t) := (f>c,tT^o4>r,-tU G H^{C{t)) where 7?.o: —?> H^{Co) is the right 

continuous inverse of the trace operator. Note that 

ll^llffi(c(t)) — ^0 \\T^o4>r,-tu\\Hi(Co) — < C 2 ||^r,-tM||ppi/ 2 , 2 (r|,) < C 3 ||M||ppi/ 2 , 2 (rp)). 

Also, we have TtU = Tt4>c,t'Tio4’r,-tU = (pr,tToT^o4’T,-tU = u by Lemma 3.5. So the function 77 = £tu — U G 
H^{C{t)) can be taken as an admissible test function in the weak formulation for £tu, and doing so yields 

[ \'^ g{t)^tu\ = j g{t)^tU^g(t)U < — / |Vg(t)£tM| H — ||rt||ppi/ 2 , 2 (p) . 

JCit) Jc{t) ^ JC{t) ^ 


□ 


Lemma 3.9. There exist constants Ci, C 2 > 0 such that for all t G [0,T] and all u G H^^^(r(t)), 


||^3(i)'^-R'‘^llL2(c^(t)) < C*! Ik “lliyi/ 2 . 2 (p(t) + ^2 Ik ^IIl2(p( 4)) + ' |r|. 

Proof. Suppose m = 0 and let 77 = £ii,tu— ^^£tu G P[l{CR{t)) which we take as the test function in the weak 
formulation for £R^tU'. 

f g(t)£R.tu\'^ < I [ \Vg^t)£R,tu\‘^ + ^ git) {{R - y)£tu)\^ , 

JCRit) ^ JcR(t) -K 

and this gives 

/ |Vg(t)£’_R,tM|^ — ^ / 4|Vg(t)£t7r|^ -I- < Cl lkllivi/ 2,2 + lkllL2(r(t)) > 

JCRit) JcrU) R K ^ 

where we replaced the integral over CR(t) by one over Clt) (this is why we need u = 0) and used Lemma 3.8 
and (11) in conjunction with (Aa). The case follows from the above and 


It7|2 


/ git)£ R,tu\^ = 

JCRit) JCrU) 


^git)£R,t{u -u) + ^Vg(t)(i? - y)u 


3.2 Truncations 

Let r be a smooth hypersurface. Define the truncation : K. - 
Tk{x) = m.ax{x, —k) + min(x, k) — x = 


<■ R at height k by 
A:sign(a;) : |a;| > fc 

X : kl < k. 
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Both Tk'. i^(r) —?> and T^: -ff^(r) —?> -ff^(r) are bounded continuous maps. The Lipschitz nature 

of the function max(-,0) implies that —>■ is bounded and by Lemma 3.3 it is also 

continuous. Furthermore, the chain rule for weakly differentiable functions u gives 

^iTku{z)) = X{\u{z)\<k}{z)^u{z) 

for almost every z. See [22, Lemma 2.89]) and the discussion after Theorem 4.3.6 in [21] for these facts on a 
domain 17. 

Now we discuss truncations over cylinders. Suppose / G C'^(K) with /' bounded and /(O) = 0. The 
chain rule Vg/(z;) = f'{v)VgV for v G H^(C) can be proved by the standard argument: approximate v by 
Vn G 27([0, oo); 27(r)), prove the identity for and pass to the limit using continuity of /' and the dominated 
convergence theorem (DCT). This then allows us to show that 

VgU+ = X{„>o}Vg'u 

(almost everywhere) by approximating r i—(r)’*' by functions with bounded derivatives, the chain rule 
and then the passage to the limit in the approximations (see [35, Lemma 1.19]). This will imply that if v and 
w are in H^{C), then max(u,w) G H^{C) and 


Vg max(u, w) 


VgU : ii V > w 
'S/gw : otherwise. 


Since u = u’*' — we have Vuj{„=o} = 0 almost everywhere. Also, if Vn,Wn are such that —>■ u and 
Wn —>■ w in H^{C), then iLnax{vn,Wn) max(u,u>) in H^{C) [35, Lemma 1.22]. Therefore, T}.{v) G H^{C) 
whenever v G H^(C). Furthermore, Tk{v) — >• u in H^{C) as k ^ oo and Tfc: H^{C) H^{C) is continuous. 

If u G then Tk{v) G Lj^i^c) 4>C-tTkiy{t)) = Tk{4>c-tv{t)) and Tk{v) G L^(0,T; iJ^(Co)) 

whenever v G L^(0, T; iL^(Co)). 

Clearly, all of this applies if we replace C with Cr and in that case we can drop the requirement /(O) = 0. 


3.3 Integration by parts 

We will need the following integration by parts result which is comparable to a result in [32] and [29, Lemma 

7.1]. 

Lemma 3.10. Let /: R — >■ K with /(O) = 0 satisfy either 

(A) / is and Lipschitz, or 

(B) / is C° and piecewise with /' = 0 outside a compact set AT CC K, 

and define F{s) = fg f(r) dr. Then for all u G W(IF^/^’^, the following formula holds: 

[ FiuiT)) - [ F{uo) - f 
/r(T) JVo Jq JT(t) 


f {u{t),f{u{t)))=[ F{u{T))- ( F{uo)- [ [ F(u(f))Vr(t) • w(f). 

Jo Jr{T) Jfo Jo Jr{t) 


(24) 


Proof. We begin with case (A). Let Un G W(IF^/^’^, L^) fl be such that Un ^ u in W(VF^/^’^, W 
Note that F(un) G L^). To see this, observe that 


\F{s) - F{t)\ < l/(r)l < ll/'ll^ maxds], JtDls - t\, 


so for almost all t, lF(it„(7))l^ < JF] 1|/'1|^ l|wn(7)|lt=o(r(t)) and since ||un(i)llL~(r(t)) bounded almost 
everywhere by 1 |m„1|^oo , we have F{un) G L^ 2 . We also see that 


\F{Un{t,x)) - F{Unit,y))\ < ll/'l|oomax(lu„(7,x)l, \Un{t,y)\)\Un{t,x) - Un{t,y)\ 


(25) 
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which shows that F(un) G ^^ 1 / 2,2 since Un € Likewise, d*{F(un)) = f{un)un G L| 2 - This means 

that the transport theorem is valid and the desired formula (24) holds for the and now we must pass 
to the limit in n. For almost every t, for a subsequence, M„(t) —>■ u{t) in LF^/^’^(r(t)), so by Lemma 3.3, 

\\f{un{t)) - /(M(t))||^i/2,2(r(t)) ^ 0 and 

\\f{un{t)) - /(w(t))||^i/ 2 . 2 (r(()) < 2 ll/'ll^ (^\\u n(^)llwi/2,2(r(t)) + lk(^)llrvi/2.2(r(t))) • 

The right hand side converges to 4 ||/'||^ ll'“(^)ll^i/ 2 . 2 (r(()) whilst the integral of the right hand side converges 

tom'wiioMmW, 2 (r(t)) since t u in L^i/ 2 , 2 - Then the generalised DCT gives f{un) —>■ f{u) 
in L^i/ 2 , 2 - For the remaining terms, we can use (25) in conjuncation with max(|a|, |5|) < |o| + |6| and 
C auchy-S chwarz. 

For the case (B), note that f{u) G L^i/ 2,2 whenever u G L^i/ 2,2 as / is Lipschitz (since /' is bounded 
a.e. and / is absolutely continuous). Given u G W(fF^/^’^, there exist Un G L^) such 

that —>■ u in W(LF^/^’^, We have that F{un) G W(LF^/^’^, L^) because F is Lipschitz. So then 

we can use the standard integration by parts formula to obtain the desired formula for Then again we 
need to pass to the limit. We have that t n in C°i by W(VF^/^’^, LF“^/^’^) ^ C °2 G^i. This implies 
that F{un) —>■ F{u) in C°i because F is Lipschitz; this takes care of the right hand side of the formula. To 
finish, since / also is Lipschitz, f(un) —>■ f(u) in L^i/ 2,2 due to the same reasoning as before. □ 


4 The harmonic extension problems on evolving spaces 

In this section, we shall consider (7) and also the following truncated harmonic extension problem”: 

given u G L^i/ 2 , 2 ) dud v G such that 

AgV = 0, TR,y=ov = u, Tfi^y^RV = 0. (26) 

As explained in the introduction, we study these problems in order to derive measurability in time of £t and 
£R,t which we do not automatically get since each £t and £R^t is constructed individually in time. 


4.1 The harmonic extension of m G L‘^^/ 2,2 

Lemma 4.1. For every u G L^i/ 2,2 with u{t) = 0 for a.a. t, there exists a U € with TU = u 

and U{t,y) = 0 a.e. t and for all y. 


Proof. Define 


U{t) 


4>c,t 



( </>r,-tM(t) \\ 

V )) 


which satisfies U G L‘j^i(^c) since (j)c,-tU{t) G L^(0,T; 7J^(Co)) by smoothness on Jq and by using (11) and 
(12) (measurability can be inferred from considerations of Nemytskii maps [33, §3.4]). It is easy to check that 
U verifies the desired properties using Lemma 3.5. □ 


Theorem 4.2 (The harmonic extension problem in the space L^). There exists a map E: L^i/ 2,2 —t L^xiC) 
such that given u G L^i/ 2 , 2 ) "v = En is the unique weak solution of (7) satisfying Tu = u in L^i/ 2 , 2 j 



Vg(t)n(t)Vg(t)r 7 (t) = 0 for all rj G L\ji{C) with Trj = 0, 


(27) 


and /p^j^(Eu)(t) = u{t). When u{t) = 0 for a.a. t, we write the solution as En. The map E satisfies 
Eu = E(m — u) + u. 
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Proof. First, suppose that u{t) = 0 for a.e. t. Let us transform the equation to one with zero initial trace. 
By the previous lemma, there exists aU € ^^i(c) with TC7 = u and crucially U{t,y) = 0 for a.a. t and all y. 
Set d := V — U G which satisfies 


Agd = —AgU and Td = 0. 


The space X := {d £ | Td = 0 and d{t,y) = 0 for all y and a.a. t}, being a closed linear subspace of 


(thanks to the continuity of T and y i—>■ d(t, y)), is a separable Hilbert space. Define J: 


: X 


by 


J{d) = i / / |V,p)d(t)|2 + / / V,q)C/(f)V,q)d(t), 

^ do dc(t) do dc(t) 

and observe that J is coercive through the use of Poincare’s and Young’s inequalities. Since J is also contin¬ 
uous, by [28, Theorem 5.25], J has a unique minimiser d satisfying J'{d,w) = {Vgd + XgU,Vg'w)i ^2 = 0 

(C) 

for all w £ X. Recalling u = d -|- Cd, we find that v £ "'^hh Tu = u and w = 0 satisfies 




{t)Xg(^t)w{t) = 0 for all w £ X. 


Jo Jc(t) 

To remove the mean value condition on the test functions, let y £ with Trj = 0 and test with 

w{t) := r]{t) — r]{t) (this satisfies Tw = 0 and w{t) = r]{t) — r]{t) = 0, so is admissible): 


0 = 


/o dc(t) 

f 

to Jc{t) 


Vg(t)U(t)Vg(t)?7(f) - 


/o ^0 


|r(t)| 


/ ^(t)) dy( [ y{t)] 

dr(t) J \Jr{t) J 


V3(i)U(t)Vg(t)?7(t) 


since v{t) = 0 for a.a. t and all y. This settles the problem for the case u{t) — 0. For general u £ T^i/ 2 , 2 ) 
define Eu := E(it — u) + u£ L‘^x{c) '''^bich satisfies jrlqj /pq) Eit(f) = u{t) and 


/o dc(t) 


Vg(t)(Eu)(t)Vg(t)y(t) = 


0 Jc(t) 


Vg(t)(E(u — u))(t)X = 0 for all 77 £ L^h^{c) T 77 = 0. 


□ 


We need to elucidate the link between E and the family of maps {£t}t^[o,T\ from Definition 3.4. 

Lemma 4.3. Let u £ L^i/ 2 , 2 - For almost all t, (Eu)(t) = £tu{t) in X{C{t)). 

Proof. Pick if £ C^{0,T) and uq G id^(Co) with TqVq = 0, then if(j)c,tVo £ ''^’bh T{ip(j)c,tVo) = 0, so it 

is an admissible test function in (27) and testing with it gives 

/ Vg(t)(EM)(f)Vg(t)^c tVo = 0 for all vq £ H^{Cq) with TqVq = 0 , for almost all t. 

Jc{t) 

By the homeomorphism properties of 4>c,t, this is same as 




for all Vt £ Pl^{C{t)) with TtVt = 0, for almost all f. 


JC(t) 

and since also Tt(Eit(t)) = u{t), we have (Eu){t) = £tu{t) by the uniqueness in Theorem 2.5. □ 

Thanks to the the previous lemma, we can use the bound (11) and Lemma 3.8 in conjunction with the 
eigenvalue estimate {A\) to obtain the next result. 

Corollary 4.4. We have for u £ L^i/ 2,2 with u = 0, 

I|E«IIl^ <C\\uh. and ||V,Eu ||^2 <C\\u\\^. 
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4.2 The truncated harmonic extension of m G fv^i /2 2 

Theorem 4.5 (The truncated harmonic extension problem in the space There exists a map 

E/j: T^i/ 2,2 that given u G T^i/ 2 , 2 ) E/jit is the unique weak solution of (26) satisfying 

Tfi,y=or’ = u and fR,y=R = 0 in T^i/ 2,2 and 

[ [ Vg(t)u(t)Vg(t) 77 (t) = 0 for all p € ( 28 ) 

Jo JCnit) ° 

and jp^Tyj /r(t)= u{t). When u{t) = 0 for a.a. t, we write the solution as ¥.ru. 

Proof. We transform (26) to having zero boundary conditions by setting w = v — {R — y)'Eu/R G 
then 

Ag(t)w(t) = - — Ag^t) {{R - y)^u{t)) on Cfi(t) 
w{t, X, 0) = w{t, X, R) = 0, 

which, by Lax-Milgram, has a unique solution w G n satisfying 

^o\yR) 

[ f Vg(t)n;(t)Vg(t)7?(t) = - / / Vg(t) [:^^ft(w(t))') Vg(t)7?(t) Vr? G 

Jo JcrU) Jo JcrP) \ R / ’ 

Indeed, define the bounded and coercive bilinear form a: L^„i , x , —>• R by the left hand side of the 

above equality and define l\ ®. by the right hand side, which is a bounded linear functional due 

to (11) and Lemma 3.8. It follows that Eru := v := w + {R — y)E{u)/R G satisfies TR^y^ov = u, 

^R,V=R = 0 and (28). □ 

Lemma 4.6. Let u G For almost all t, {ERu){f) = SR^tuif) in H^{CR{t)). 

This lemma follows just like Lemma 4.3 since (fc.t ■ Hq{Cr{0)) —>■ Hl{CR{t)) is a homeomorphism, and it, 
along with Lemma 3.9, implies the following. 

Corollary 4.7. There exist constants C\, C 2 > 0 independent of R such that 

||E_Rn||^2 < Cl 11^11^2 + 2\JR\V\ ||u|| 

l'YCr) 

||VgEfln|| 2 <C' 2 ||m|Il 2 */i?>l. 

" '^wi/ 2.2 

A third way to interpret the map !^r from (18) is as a map ^r \ {77 G | TR^y^Ri] = 0} —^ R^h^{c)’ 

and again this is an isometry. 

Lemma 4.8. We have SCrEru —>■ Eit in 

Proof. Lemma 2.13 gives for almost all t \\yg(t) 1 ^r£ R,tu{t) — Vg(t)ftu(f)||^ 2 (c((p 0 and to use the DCT it 

suffices to find an integrable uniform in R bound on the above norm which follows from Lemma 3.9. □ 

5 The non-degenerate problem: proof of Theorem 1.4 

Let /3: R. —>■ R be a function satisfying {Ap) on p. 5. We will prove Theorem 1.4 in this section, that of the well- 
posedness of problem (P/3). For easier reading, we will shorten the duality products (•, •)w-i/ 2 , 2 (rp)) jyi/ 2 , 2 (r(()) 
to (•,•) (an abuse of notation) and (•, •)^^-i/ 2 , 2 (rg) vvi/ 2 . 2 (r(,) to (•,-)o. 
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5.1 Existence of solutions to the truncated problem 

This subsection is devoted to the proof of the following theorem. 

Theorem 5.1. For each R>1, there exists a unique weak solution uh G W(1T^/^’^, to (P/ 3 fl) with 

VgEfl(/3(ufi)) G and ur{0) = uq satisfying 


/ {'aRit),v{t)) + / Mfl(t)r?(t)Vr(t) • w(t) + / / V-g^t)8H,t{P{uR(t)))Vg^t){ER'q){t) = Q 

/o do JT(t) do dCB(t) 


for all rj G T^i/ 2 . 2 , w/iere Ert] G satisfies TR^y^oERiq = g and TR^y=RERr] = 0. 


We hide the subscript R in ur and write just u for simpler notation. Define a_R(f;-,-): x 

W^/^^‘^{r{t)) ^ R by 


aR{t;u,r]) = 


ICrU) 


^ g{t)^ R,til3iu))V g(^t)ER{t)r) 


where ER^t): —>• E[^{CR{t)) is an (arbitrary) extension that satisfies Tr t y=oiER{t)ri) = rj and 

TR,t,y =R{ER{t)vi) = 0; the choice of Er does not matter (see Remark 2.12). To derive the Galerkin approxi¬ 
mation, we pull back the first two terms in the equation onto T 0 and then make a substitution in order to put 
the Jacobian term onto the elliptic form. Let {bj} be an orthonormal basis of L^(ro) that is orthogonal 
in lT^/^’^(ro) and let u{t) = (j)r.-tu{t). The Galerkin approximation is the system 


/ w'nW&i + / Un{t)bj(j>_t{^Y'(t)-^{t))+aR{t;(j)tUn{t),Jo(j>tbj) = 0 Vj = l,...,n 
Jto Jto (29) 

M„(0) = Uon 

for an ansatz Un{t) = oti{t)bi with unknown coefficients Ui = of and uon G Ei(0) := span{6i, is 

such that uon -t wq in and ||uon||wi/ 2 , 2 (ro) < C\\uo\\^i, 2 , 2 (Tg)- 

Remark 5.2. We pulled back the equation onto a reference domain in order to facilitate the procurement of a 
bound on u'^ which is needed for a strong convergence result. This transformation to the reference domain Tq 
could have been avoided if we knew that the orthogonal projection operator Pf: L^(T(t)) —>■ V„{t) := 0t(Ei(O)) 
defined by 

{Pfu - u, fn)L2(r(t)) = 0 far all Vn G 14(t) 

is bounded as a map Pf:V{t)^V(t). Such a bound is true when t = 0 because of the special choice of 
basis functions, but for arbitrary t the desired bound appears elusive. Of course, such a result would be of 
fundamental use generally in evolutionary equations on evolving domains. 

Lemma 5.3. The Galerkin equation (29) has a solution Un G i?^(0,T; td„(0)). 

Proof. The equation (29) leads to 


Up n 

ai{t) / 6 i 6 j(/)_t(Vr(i) • w(t))-k aK(t; 

i—1 i—1 

aj(0) = (uon, 


Vj = l,...,n 
Vj = l,...,n. 


(30) 


Define a.(t) = (a^t),...,a„(t))T, b(t) = (^* 61 ,..., (/)t 6 „)T, a(t, a) = {aR{t-,a ■ b{t), ... ,aR(t-,a ■ 

b[t), J^cftbufy and the matrix (W{t))ij = Jj,^bjbi(j)-t{\7r{t) ■ w(t)). The system of equations (30) is then 
written 


a'{t) = F{t, a{t)) := —W {t)a(t) — a(t, a(t)) 
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with initial data Q!(0) = ((uon, &i)L2(ro) ; (uon, bn)L^(ro))''^- We need to show that t ^ F{t, a) is measurable 
for fixed a G R". The term with the matrix is clear. For the remaining term, we have 


aR{t;a- b(t), Jl(j)tbj) = 


lo Jr(t) 

rR 


y -g(t)£ R,t{P{oL ■ b{t)))ygi^t)ER(t){Jl(j)tb:i) 




^0 Fq 

+ r [ J^tdy^-t£R,mo,-b{t)))dy<j^.tER{t){4<j>th,), 
Jo Jvq 


and we know that £r4P{oc ■ b{t))) = EH(,5(a • b)){t) for almost all t (Lemma 4.6), and the pullback of the 
latter is measurable as a function of t since E/i(/3(Q: • b)) G the same argument can be used to deal 

with the ER{t) term. Now suppose that —>■ a in R". We see that 


i(t,a^) - a(t,a)||^„ = ^ \aR{t]a^ ■ b{t), J^cjjtb,) - aR{t;a- b{t), jQ(j)A)\''" 




lo Jr(t) 


■ b{t)) - /3(a ■ M^))]Vg(t)£^/j(t)(To</'t&*)l 


< C{R) \\P{a^ ■ b{t)) - I3{a ■ 6(t))||^i/2.2(r(t)) E 


by Lemma 3.9, and this tends to zero by Lemma 3.3 since a.^ ■ b{t) —>■ a • b{t) in Therefore, 

a i—^ a(f, a) is continuous and so F is a Caratheodory function. The uniform bound that we shall derive in 
the next subsection shows that ||Q:(f)||jj„ < c for all t if ct satisfies the ODE (30). Let us now prove that there 
exists / S L^(0,T) with ||F(t, a)||jj„ < f(t) for every a S {a G R” | ||Q;|lRn < 2c}. Note that 

*^)IIk" — E/ \\'^g{t)£R,til3io^ ■ b{t)))\\j^2(^CR{t)) \\'^9it)^R(04o4>tbj)\\j^2(^CR{t)) 

3 

< C*! E ■ ^W)ll^i/ 2 , 2 (r(t)) \\'^-g{t)ER{t){Jl(j)tbj)\\\ 2 ^CR(t)) Lemma 3.9) 

3 

< Cl ii/3'iiL ii^iir E wumwy^^Hnt)) E l|VM*)C«(0(4</>t^i)|ll2(c,p)) 

^ 3 

< C 2 ||ci|lKn ll^t^i|lwi/2.2(r(t)) \\'^git)ER(t){jQ(j)tbj)\\j^2(^CR{t)) 

i,3 


SO that overall (because the Frobenius norm H-Hp is compatible with the Euclidean vector norm), 


||C(t,a)||K„ < 2c \\W{t)\\p + vF^E ll'?^i^*llivi/ 2 , 2 (r(t)) ||^§(t)CflW(-/o<(>i&i)|L 2 (c„p)) =: fiO 


\ bi / 

and the term in the brackets on the right hand side is integrable over (0, T). Now an application of the ODE 
theory in [62, Problem 30.2] gives global existence of a solution u„ : [0, T] —>• 14(0). □ 


5.1.1 Uniform estimates (in n) 

Multiply the first equality in (29) by aj{t) and sum up to get (using an arbitrary linear extension) 


Now, in 



U„{t)Un{t) 


■ w(t)) + aR{t] (j)tUn{t), Jo4’tUn{t)) = 0. 


^r{I'-! ^tRn{i): Jo^tRnib)) — / ^ R,t{P{4^tRn(l)))^ g[t)ER{t)(^jQ(l)t'Ujn{t)) 

JCait) 
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let us pick Eii{t){Jl(j)tUn{t)) = Jo/3 ^{£R,t{l3(4>tUn{t)))), which is valid since TR,t,y=oER{t){Jl(j)tUn{t)) = 
Jll3~^{P{(t>tUn{t))) = Jo(j)tUn{t) and TR,t,v=RER{t){J^(j)tUn{t)) = 0, and this gives 


1 d 

2 dt 


To 




< 2 ll^r • w|l^ 


wherein we note that 

rR r 


'To 


JCait) 
i(t)^ + C 2 


\y g(t)£ RAPi^tUnit)))^^ 


ICR{t) 


C,\/ 3 -A£R.tWtUnm)\ + e\Vr£R,tWtUn{t)W 


/ / \rA£R,tmtunmr < cs iKr^/L mAunirnihim) < ^4 wnAmhivo) 

»/ 0 

by (14), where C 3 and C 4 depend on R and Ai. Then Gronwall’s inequality implies 

inax ||un(t)|li 2 (r„) + ||VgE/j(^((()(.)U„)) ||^2 < C. 

te[o,Tj l^{Cr) 

Remark 5.4. ITe needed to truncate the domain in order to obtain the previous bounds. If the domain was 
instead the full cylinder C(t), the extension of the test function would have to include a cut-off function so 
that it belongs to for example, ififp is as in Definition 2.8, then we could choose 


E{t){jQ<ptUnit)) = Jo/3 ^[£t{l3{(l)tUn{t)) - AcftUnit))) + f)p^cjitUnit))] 
but this leads to a residual term of the type 


P ^AtiPiPtUnit)) - PiptUnit))) +lljpPi(l)tUn(t))]Vr£tiPi(l>tUn{t)))VrJo 

Jo Jr(t) 

and we would have to make restrictive assumptions on the evolution to neglect this term as we send p —>■ 00 . 
Writing P{un) = Tfl:,j/=oEfl(/3(ura)) and using the trace inequality, 


||/3(u„)||i2 < Cl ||Efl(/3(u„)) 

Wl/2,2 " 


1^2, <C2||/3'C|K||i2^ +C3<C4 


'' <C2\\P'f 

"h^{Cr) 

by Corollary 4.7 (C 2 and C 4 will depend on R) and the energy estimates. Since P~^ is Lipschitz, this implies 


llWnllj 


Vi/ 2,2 


< c 


independent of n (using the boundedness result of Lemma 3.3). The bound on the time derivative follows 
too: take 77 S lT^/^’^(ro), recall (29) and that is self-adjoint: 

/ Ru{t)v = - [ u„(t)P°(p)(/)_t(Vrp) • w(t)) - / f Vg^t)£RAPiPtUu{t)))Vg^t)ER{t){4(j)tPAv)) 

JTn Jo Jr(t) 


assuming a linear extension. Now picking Eii{t)(jQ(j)tPn{v)) = -^o^R,t{4>tPni'n))i observe that 


\Vg^A4£RAAPum\\R.^eRit))<c;4 


lo Jrit) 


\Wg^t)£RAPtPAv))f + \£RAPtPAv)A 


< C2II(^) II 14/1/2,2 

(where again C 2 depends on R and Ai) which implies that 


(by Lemma 3.9 and (14)) 


/ {Rn{^)y V)o — C 3 II P„ II L2(0,T;W'-/2>2(ro)) ( L'^{ 0 ,T-,L'^{ro)) || ^ff(t)®^fl(/3(<?^t'art)) || p2 

^ C '4 ||?7||p2(o_T;rvi/2,2(ro)) 

by using the uniform estimates. Now taking the supremum over rj € L^(0,T; tT^/^’^(ro)) shows that 


l^nllL 2 (o^T;ll/-i/ 2 , 2 (Pp)) < C. 
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5.1.2 Passage to the limit in the Galerkin approximation 


We have as n —>■ oo 


Un^U 

in L2(0,T;lU^/2’^(ro)) 

U'u u' 

in L^(0,T; 

Un ^ U 

in R{Q,T-L\To)) 

D^Rli{Un) 0i 

in ^‘l'^(Cr) 

dyERP{Un) 0y 

in Rl'^(Cr) 


( 31 ) 


where D_^ = (Vr)i is the z-th component of the tangential gradient and Aubin-Lions yielded the strong 
convergence. Therefore we have —>■ m in L ^2 and P{un) —>■ /3('a) in thanks to the Lipschitz continuity 
of /3. Using the boundedness of E/j in the space L^i^ 2 ^Cr) Corollary 4.7, we obtain by linearity that 


Ei{(/3(M)) in L 12 (^Cr)- 


(32) 


Identification of the spatial term Take the test function 

r]{t,y,x) = tp{t){(j)tVo){x)h{y) where z/> G C))°(0,T), vq G ^^(ro) and h G C^{0,R) (33) 

belonging to in the spatial integration by parts formula [ i, §2.1] on r(t) integrated over y and t: 


10 JCR{t) 


{DiEji/3{un))r] = - 


10 JCrU) 


(EK^(M„))i2i?? 


/O JCR{t) 


{Eiil3{un))vHi'f 


where H is the mean curvature. Using (31) and (32), we have 


/o JCR{t) 


OiV = - 


0 JCR(t) 


{Eiil3{u))D^ 


lo JcrU) 


{ERji{u))r]Hv^, 


whence it follows that for almost every t, for almost every y, 


/ Si{t,y)(j)tVo = - (Eiil3{u)){t,y)D^(j)tVo 

lr{t) Jr{t) 


'm 


{ERP{u)){t,y)(j)tVoH{t)vf {t). 


Since this holds for all <j)tVo G (^'^(ro), it also holds for all v G C'(](r(t)), which implies that Di{ER/3{u)) = 9i 
by definition. 


Identification of the y term Again take rj G L,^hi(^Cr) consider the integration by parts 

formula 


lo JCR{t) 


{dyERl3{un))y = - 


lo JCR{t) 


{ERl3{un))dy'q. 


As before, passing to the limit we find for almost every t and almost every x (again since (jjtVo ranges over all 

of Ci(r(t))) that 

pR pR 

/ 0y{t,x)h = - (ERP{u)){t,x)dyh 


and thus we identify 0 y = dyERfii(u). 


Conclusion of the proof of Theorem 5.1 Therefore, the last two convergences listed in (31) can be 
replaced with VgERl3{un) VgER/3{u) in L^ 2 (Cjj)- Recall that Vn{t) := span{(()t&i,..., (/)(&„}. Given y G 

L‘^ 1 / 2 , 2 , by density, there is a sequence r]i{t) = X)j=i lj{t) 4 >tbj with ry G Ly^ such that ry ^ y in L^i/2,2- If 
I < n, then yi G Ly^ and we multiply (29) by 7j(t) and sum up to get 



u'„it)m(t) + 



Mn(t)?7i(t)<(’-t(Vr(i) • w(t)) + 



Vg(t)£:i^,i(/3(■u„(<)))Vg(t)Ui^(^)(J*?7^(^)) = 0, 
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where fjiit) := We obtain after integrating the above equation and sending n —>■ oo the equation 




' 0 
pT 


u(t)i?i(t)0_t(Vr(t) • w(t)) 


+ 


/O JCnit) 


Vg(t)£RAP{u{t)))yg(t)ER{t){Jlr^i{t)) = 0. 


(34) 


Let us prove that (j)-t{Jo)vi 4>-tiJo)fi in L^(0,T; W^/^’^(ro)). For the seminorm, we have 


r 


\(l>-tiJo)ix) {[fjiit) - fi{t)]{x) - [fjiit) - fjit)]iy)) 




'-y|' 


iTo JTo 


-y\^ 


< Cl - r7(t)|^i/2,2(ro)+(^ 2 ^ \[fjiit) - fjit)]iy)\^ j^^—^^^dxdy 


{(p-tJlj is Lipschitz with uniform Lipschitz constant) 


< Ca Wmit) - ?7(t)Ci/2,2(ro) 


see Lemma 5.2 in [4] for last line. Integrating over time and passing to the limit shows the result. Thus 
JoVi —t JqV in L^i/ 2,2 and it follows from Corollary 4.7 that VgEfl(? 7 iJg) ^ VgE,iiir]jQ) in With 

this in mind, taking limits / —> oo in (34) with En = E/{, and then because the extension can be arbitrary 
(Remark 2.12), we obtain 


{u'it),fjit))o+ [ [ M(t)r 7 (t)(()_i(Vr(t) • w(t)) + / f \/gy)£R^til3iuit)))\7gy)ERiJ^'\)it) = 0 

Jo Jfo Jo JcR{t) 


for all fj € L^(0,T; IF^/^’^(ro)). Now, pushing forward the first two integrals, recalling from the proof of 
Theorem 2.33 of [3] that iiit) = (/)1(( J°i2'(t)) (where is the adjoint of </)_(: IF^/^’^(r(t)) —>• IT^/^’^(ro)) 
and using 


{fi'it),f,it))o = {j°u'it),iJ^)-^fiit))o = {r-tiJ?u'it)),MiJ?r^m)) = (n(t), j^</>t^(t)) 


gives 

f mt),J*(l)tfiit)) + [ [ uit)(j)tfjit)Jl\7r(t)-^{t) 

Jo Jo Jr{t) 

+ r [ Vgy)£R,muit)))S/gy)ERi4U(-)m) = 0. 

Jo JcrU) 

Picking fjit) = (j)-tr]it)/(j)-tJQ yields that u satisfies the equality given in Theorem 5.1 for each y G L^i/ 2 , 2 - 
A standard argument involving integration by parts of the equation satisfied by u and the equation satisfied 
by Un and then passage to the limit in n shows that the initial condition is satisfied, see [3, §5.3]. 

5.2 Existence of solutions to the non-degenerate problem 

Therefore, for each i? > 1, we have a function ur G W( 1 F^C> 2 ^ ^- 1 / 2 , 2 ^ m_r( 0) = ug and VgEi{(/l(ufi)) G 

^h(CR) satisfying 


{uRit),yit)) + [ Mfl(t)??(t)Vr(t) • w(t) + / Vg(t)fii,t(/3('Ufi(t)))Vg(t)(Aii7?)(t) = 0 (35) 

Jr{t) JcRit) 

for all y G L'^-^^ 2,2 and for almost all t. We now want some estimates independent of R. 
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5.2.1 Uniform estimates (in R) 

Let us pick r] = uji and use Ert] = /3“^(Efl(/3(u/j))): 

— [ J • w(t) + r [ (/3-i)'[ffi,,(/3(ufi(t)))]|V,p)ffi,,(/3(ufi(t)))|2 = 0. 

^ Jr{t) ^ Jr{t) Jo Jr{t) 

Since {/3~^y > Cpi_ ^ and ur{0) = uq, we immediately obtain via Gronwall’s inequality that 

max ||nfl(t)||i2(r(t)) + ||Vglfl(/3(ui{))||^2 < C 

tG[0,T] L^(Cji) 


independent of R. Choosing Er = 
supremums over 77 G L^i/ 2,2 gives 


E/{ in (35), using Corollary 4.7 and the uniform estimates, and taking 




L 


2 

lv- 1 / 2,2 


< c. 


Lemma 2.10 implies that |7tw|i,j/i/2,2(p(j)) < C ||Vg(t)'t'||p 2 (c(tp all v G we claim the constant is 

independent of t. Indeed, an inspection of the proof of the lemma reveals that we need to check whether the 
trace map 7): H^{C(t)) —)> W'^/^’^(r(t)) is bounded uniformly in t and whether the constant in the Poincare 
inequality on r(t) is independent of t. The first question has been settled in §3 and the second is also 
affirmative due to [4, Lemma 5.9]. Using this inequality, we find 




which implies that \/3(ur)\]^2 < C. This gives boundedness of ur in the fractional seminorm, and thus 




^ 1 / 2,2 


< C. 


5.2.2 Passage to the limit in R 

Therefore, we have 


UR 


u 

in 

T ^ 

^VVi/2,2 

UR 


it 

in 

^10 

1 

to 

to 

UR 


u 

in 

IjL ‘2 

i2j^flEfl(/3(u_R)) 


9^ 

in 

T ^ 

^L-^{C) 

dySfRERifiiUR)) 


9y 

in 

T ^ 

^L-^{C) 


and we need to identify the limits 9i and 9y. Our first task is to show that 3fRER[(i{uR) — I3{ur)) — >■ 
E(/3(u) — /3(u)) in L‘^ 2 (^Qy Set wr = P{ur) and w = P{u); since WR{t) — WR{t) —>■ w(t) — w(t) in L^(r(t)) for 
almost every t, by Lemma 2.15, 


fR{t) ■■= 


3fR£R{wR{t) - WR{t)) - £{w{t) 


w{t)) 


2 

0 . 

L 2 (c(t)) 


By virture of ^ being an isometry, and using (11) and (14) 

2 

LHm) 


|/ij(i)| < Cl WR{t) - WR{t) 


w(t) — w{t) 


L 2 (r(t)) 


=: 9R{t). 


LRr(t)) 


Now, for almost all t, gR{t) —>■ 2Ci w{t) — w(t) 
generalised DCT, lim/{_>.oo Jq fR(t) = /o^linifl->oo fR(t) = 0 , giving 

^rEr(wr - W^) E(w - w) in Lp 2 (c) 


and /o 5fl(t) 2Ci\\ w — rc||p 2 , so that by the 


1,2 
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as desired. Now, choosing rj as in (33) except with h G we have 

f [ Di^nEHWRT]=f f Di^R (eh{wr - W^) + ^ V 

Jo Jc(t) Jo Jc{t) \ -K / 

= - [ [ D^r]3fRER{wR-w^) + ( [ 2fRER{wR-w^)r]Hi'f 

Jo Jc{t) Jo Jc{t) 

and passing to the limit on both sides gives 

rT P pT 


[ [ = — [ [ DiV^iw — w) + [ f E{w — w)r]Hvy, 

Jo Jc(t) Jo Jc(t) Jo Jc(t) 


T 

0 Jc{t) Jo Jc(t) Jo Jc(t) 

and then an argument similar to that in §5.1.2 shows that D_jE{w — w) = 0i. For the y term, 

_^ , R-y _\ 


[ [ dy^RERWRr]= [ [ dy3fR(ER{wR-WR)-\ - ^wr \ rj 

lo Jc(t) Jo Jc{t) \ ^ J 

pT p 


_ Wk 

2 ^rEr{wr - WR)dyr] + Xv<R-^V, 


lo Jc{t) 


R 


(36) 


where the last term on the right hand side 


/ / Xy<R^V = [ / Xy<R[y)Ky) [ (t>tVo -)■ 0 

Jo Jc(t) Ji Jo ^ Jo Jr{t) 

since WR{t)'il^{t)/R —>■ 0 and Xy<R{y)^{y) /o°° ^(y) both due to the DCT (recall that WR{t) —>• w{t) 

a.e.). Then taking the limit in (36), we get 

rT r pT r 


/ / ^y'n = - / ®:(u> - w)dyr] 

Jo Je{t) Jo Jc(t) 


0 JC{t) 


which again gives 9yE(rt; — W) = 9y by similar reasoning to §5.1.2. Now, integrating (35) in time, we can pass 
to the limit by first of all taking Ert] = SfiEiy (this satisfies i?fl77|y=o = V and ERr]\y=R = 0 since R > 1). 
Replace the integral over Cfl(t) by one over C{t): 


[ {uR{t),v{t)) 
Jo 




UR{t)r]{t)\/r(t) ■ w(t) + 




Vgp) J"fl£:fl.t(/3(nfl(<)))Vg(t) J"iEi7?(t) = 0, 


and then using the above convergence results and recalling that the elliptic form can have an arbitary extension, 
we find exactly the weak formulation (9) of Theorem 1.4. For the conservation of mass, note that Ufi(t) = 
Jp^ uq holds simply by testing with rj = Ert] = 1 and then we can use the strong convergence of ur to u in 
and the continuity of t !->■ to get the result for all t. 


5.3 Contraction principle 

Let uoi, uo 2 G L°“(ro) be initial data and consider the respective solutions uir and U 2 r to the truncated 
problem (9). The contractivity can be proved with a sensible choice of test function (for example, see [11] for 
a continuous dependence argument). Take the difference of the two weak formulations, set ViR — E/{(/3(ui fi))> 
pick 7] = ^T^{{uir — U 2 r)'^) and integrate over time: 

- [ {d'{uiR — U2R),Te{uiR — U2 r)~'') + - [ [ (uifl — M 2 fl)F'e(ui/J — M 2 fl)'^Vr(t) ' w(t) 

e Jo e 7o Jr{s) 

+ / / '^g{s)iviRis) -V2Ris))X/g(^s)ER{s)vis) =0. 

Jo J Ch( s) 
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Defining S^{s) := Tc(r~^)/£ dr, applying Lemma 3.10, taking the limit inferior and using Se(-) —>■ we 

obtain 


I (Mifi(t)-U 2 fl(t))++liminf i / / Vg,^s){viR{s) - V2 r{s))\ 7gt^s)ER{s)r]{s) 

Jr{t) <=-*-0 e Jo Jcr{s) 

= [ {uir{ 0) - U2Rm+. (37) 

drn 


Let us pick Ert] = ^(uii?)-/3 ^(u 2 ii))+)/e G which satisfies TR^y^Q{ERr]) = T^{{uir-U 2 r)"^)/ 

e and 'ii^R,y=R{ERri) = 0 so is an admissible test function. Here we used that, for example, TR,y=oT,= (w'^) = 
Te{TR,y=ow)~^ for all w G E[^{Cr{0)); this holds due to a density argument using the continuity of o (•)+ 
between H^{Cr{0)) (see §3.2) and between lL^/^’^(ro) (by Lemma 3.3). We also used that TR,y=oP~^iw) = 
(3~^{TR,y=ow) for all w G H^{Cr{0)), which follows again by Lemma 3.3 and the continuity of the map 
/3“^: E[^{Cr{0)) —^ H^{Cr{0)) (which is a consequence of the boundedness and continuity of 

Setting B^{s) := {(x, y) G r(s) x [0, i?] | 0 < /3~^{vir{s, x, y)) — /3“^(u2fl(s, x, y)) < e}, the elliptic form in 
(37) is (see §3.2) 


'Cr{s) 


^g{8){viR{s) - V2 r{s))S/ g(^s)ER{s)r]{s) 


> 


- [ ygis)(,ViR{s) - V2R{s))yg(s)W \viRis)) - ^ ^(W2fl(s))) 

« Jb,{s) 

- f ((/ 3 ”^)'(l'lfl(s)) - (/J~^)'(v 2 R(-s)))Vg(,)(viR(s) - V 2 R(s))Vg(,)V 2 R. 
e JbRs) 


(38) 


Here,/S ^(vir) -^(v 2 r) = (/3 ^)'(viR)(Vg(s)ViR - Vg(^)V 2 R) + ((/3 ^)'(vir) - (/J ^)'(v 2 R))Vg(s)V 2 R was 
used to derive (38). The right hand side of (38) can be estimated as follows: 


-/ 

e Jb^(s) 

<-||(r')'' 

e " 


((/3 )'(^^lfi(s)) - (/3 )'(l^2ij(s)))Vg(s)(uifl(s) - V2R(s))Vg(,)V2R(s) 


lB,{s) 


l^^lfl(s) - U 2 fi(s)||Vg(s)(uiK(s) - U 2 fl(s))||Vg(s)U 2 fl(s)| 


^)"lloo ll^'lloo / 1 ^ ^(?^lfl(s)) -/3 ^(v 2 R(s))\\^g(s)(viR(s) -V 2 R(s))\\Wg^s)V 2 R(s)\ 

^ Jb,{s) 


<ll(r')"ILll/ 3 'll 


'Cr{s) 


XB,is)\'^gis){viR(s) - U 2 fl;(s))||Vg(s)U 2 fl;(s)|. 


(39) 


Now we show that this expression tends to zero as e tends to zero. By DCT the integral on the right hand 
side of (39) converges to the integral of the limit, so we shall focus on the pointwise limit of Xb^(s)) namely. 


X{zGCb(s)|/3-i(i;i_r(s,z))-/3-i(«2b(s,2:))=0} X{zeCR{s)\viR{s,z)-V2R{s,z)=0} ■ 


(40) 


Observe that Vg(;,) (uifl;(s) — U 2 fl(s))|{„ij,(s)-t) 2 ii(s)=o} = 0 a.e. on [0, i?] x r(s) by a theorem of Stampacchia 
(see §3.2), so if {/3“^(uifi(s)) — fJ~^{v 2 R{s)) = 0} has positive measure, then the limit of the integral on the 
right hand side of (39) vanishes. So then let us suppose that /3 “^(ui_r) —/3“^(u2_r) = 0 only on a set of measure 
zero. In this case, (40) is exactly 0, so again the limit vanishes. This implies in (38) that 

li^nf I Vg(s)(uifl(s) - V2Ris))yg(^s)ER(s)r]is) > 0. 


'Cr{s) 


Plugging this back into (37), we obtain the desired contractivity for uir — U 2 r at each point in time. Now, 
by the work in the previous subsections, thanks to the strong L ^2 convergence, it follows that for almost all 
t, uiR(t) -X ui{t) and U 2 R{t) -X U 2 (t) in L^(r(t)) for a subsequence. Therefore we can pass to the limit and 
we will obtain for almost all t 


[ iui{t) - U2it))~^ < [ (Uoi - U02)~''■ 

Jr{t) Jfo 
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6 The fractional porous medium equation: proof of Theorem 1.2 

We pick (see [30, p. 102]) a sequence of smooth functions 'I'fe such that 4'fe(0) = 0, 'I'fe —>■ 'h in 

< Ci\r\ + C 2 , < C'fe, and 1/Cfc < < k. The previous section gives us existence and 

uniqueness of Uk G satisfying 


{uk{t),r]{t)) + [ u/c(t)77(t)Vr(t) • w(t) + / Vg^t)£t{'i’k{ukit)))Vg(^t)Eit)r]{t) = 0. (41) 

Jr{t) Jc{t) 

Now we obtain appropriate estimates independent of k and pass to the limit for the last time. We first look 
for a weak maximum principle. Let us set Wk{t) = Uk{t)e~^^ (note that Uk{t) = e^*{wk{t) + Xwk{t))) and pick 
V = {wk — M)+ where M := We would like to pick the extension of ry = {wk — M)+ to be 

(4-^i(E(4'fc(ufc) - ^;M))e-^‘ -mY 

but this is not possible since the bracketed term is not square integrable. Therefore we define 

gk{uk,p) ■= E(4'fc(ufc) - 4'fc(ufc)) + tpp'^k{uk) 

and pick E{wk — M)+ = {^^^^{gk{uk, p))e~^* — Mil>p)~^ G which satisfies Ti?77 = {wk{t) — M)+ and 




{'^k^)'(9k{uk,p))e ^^\Jg(t)gk{uk,p) - : if 4'^^(gfc(ufc,p))e ^^-Mtpp>Q 

0 : otherwise. 


Equation (41) reads 

{wk{t),{wk{t) - M)^) + [ wfe(t)(wfe(t) - M)+(A + Vr(t) • w(t)) 
Jr(t) 


/r(t) 


Jc(t) 

and the gradient term on the set {gk(uk, p))e~^^ — Mtpp > 0} is 


^g{t)£t{^k{uk{t)))yg{t)E{t){wk{t) - M)+ = 0, 


/C(t) 


^ g{t)£t{'^k{uk{t)))y g{t)E{t)g{t) 


/C{t) 


7g{t)£t{'i’kiukit))) (^{'i’f.^y{gkiuk,p))e ^*(Vg(t)E(4'fc(ufc(t)) - 4'fe(ufe(t))) + ■!/'p4'fc(wfc(l))) - 


)c{t) 


i'^k )\9k{uk,p))e [\Vg(t)^t{^k{uk{t)))y+dyEt{^k{uk{t)))'ilj'p^k 


JC{t) 

> e~^*k~^ 

1 


Mdy£t{'^k{uk{t)))ilj' 


lc(t) 


\yg(t)£t{^k{uk{t)))\ - e Ck / \dy£t{'^k{uk{t)))\\y’'p'Sk{uk{t))\ 


lc{t) 


> 7:e-^*C,ik) / \Vg^t)£ti^kiukmr - :.e-^‘C2(fc)|r(t)||vI/fc(u,(t))| 


/C(t) 


where the last term in the antepenultimate line vanished because dy£t{^k{uk(t)) — kiuk{t))) has mean value 
zero and to derive the last line we used Young’s inequality and that jV'pl < C/p with supp(p) C [p, 2p]: 


/C(t) 


\dy£t{'^k{uk{t)))\\'i/'p'^k{uk{t)) \ < 


p2p 


/C(t) 


i\Vg^t)£ti^k{ukmr+CeC\Tm^kiuk{t))y / -. 
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Thus, we have 


{wk{t),{wk{t) - M^) + ( Wfc(i)(wfc(t) - M)+(A + Vr(t) • w(i)) 

Jr{t) 

~ ^(sfc(Mfc,p))e-^‘-Mi/>p>0}® C'2(fc)|r(i)||4'fe(ufc(i))| p — 

Choosing A := ||Vr • w||^ and sending p —^ oo, we can discard the last two terms and we will find 

^4 / iiwk{t) - M)+f - i / {{wk{t) - M)+)2Vr(t) • w(<) < 0. 

^ JT{t) ^ JT(t) 

Gronwall’s inequality implies boundedness of Wk and hence 


esssup||ufe(t)||i„o(r(t)) +esssup||«'fc(ufe(t))||ioo(r(t)) < C. (42) 

tG[0.T] ^ ^ tG[0.T] ^ ' ” 


The second bound holds because esssupjgjQ-rj ll'I'fc('«fc(0)llL~(r(t)) — ^n^x (|'I'fe(C')|, |4'fe(—C)!) (since 4'^ is 
increasing) and the right hand side is bounded since ^ 4/. 

Now we focus on obtaining a bound on 4'fc(itfe) in T^i/ 2 , 2 - To this end, let us define the antiderivatives 

Hk{r) = / 4'fc(s) ds and Gk{r) = 

Jo 

and also the antiderivatives H and G by the obvious formulae. If m G L'^{M), then Gk{u) G L^{M) and 
Hk{'^k^{u)) G L^(M); this follows from 

|Gfc(M)| < max(|4'-i(u)||u|, |4'-1 (-m)|| w|) < (Ci|m| + C 2 )\u\ (43) 


(s) ds 


and 

= u'l’k'iu) - Gkiu) < \u\iGi\u\ + G 2 ) + \Gkiu)\ < (Calwl + Ci)\u\. (44) 

These properties are also true for G and H. 

Remark 6.1. ITe could have generalised the porous medium nonlinearity 4'(r) = |r|"*“^r to simply having 4^ 
as a eontinuous inereasing function. In this case 4' is no longer invertible so we would have to use Legendre 
transforms [30]. 

Test the equation (41) with rj = ^k{uk), pick E{^!k{uk)) = E(4'fc(wfe) — ^k{uk)) + ifp^k{uk) and use the 
integration by parts formula of Lemma 3.10: 


/o Jo Jvit) 


|Vg(t)ft(4'fc('Ufc) - 4'fe(u/c))C + dy8t{^k{uk) - ^k{uk)Wp^k{uk) 


< f Hk{uo) + [ [ iLfe(Mfc)Vr(t) ■ w(t) - f f Wfc4'fe(Mfc)Vr(t) • w(t) 

Jro Jo Jr{t) Jo Jrit) 


where we threw away the Hk{uk{T)) term since Hk > 0. The second term on the LHS disappears since 
the harmonic extension of a mean value zero function has mean value zero too. Then we finally get after 
using (44) that \Hk{uk)\ < Ci ||4'fc(wfe)||i~^ + G 2 . This takes care of the second term on the right hand 

side, and as for the initial data, we note that \Hk{uo)\ < Gi ||4>fe(ito)||£c=o(rg) + G 2 and ||4>fe(uo)||ic=o(rp) < 
max (|4'fc(M)|, |4'fc(—M)|), and the right hand side is bounded like before. Thus 


VgE(4'fc(ufc) - 4'fe(Mfc)) 


lHc) 


< G. 


Writing '^k{uk) = TE(4'j,(Mfe)) and using Lemma 3.6 and the previous uniform bounds, we have 


l^/c(w/c)||j 


'wi/ 2.2 


< G. 


31 













Finally, integrating and rearranging (41): 


{uk{t),T]{t)) < II Vr ' 


hk\\LK MlK + ||Vs]E(4'fe(ufc))| 


" 1.2 


"l2(C) 


INsEvW 


choosing Er] = SfpEpT] for some p > 1 and using Vg^Epp 2 = VgEpp 2 

L^iC) L'^iCp) 

the last inequality by Corollary 4.7, it easily follows that 


L'^(C) 

< CM 


^wi/ 2.2 


with 


independent of k. Therefore, we have 


\Uk\\ 


p- 1 / 2,2 


< C 


Uk ^ u 

in L|2 

Uk ^ u 

Lyy-1/2,2 

Vk ■= 4'fe(wfe) ^ V 

ill Lyyi/2,2 

VgE(nfe - Vk) ^ 9 

in -^L2(c) 


(45) 


(46) 


with the strong convergence by Aubin-Lions. Now the question is whether v = 4'(u). If so, then we can also 
identify 0: indeed, we know that the map G: i^i/ 2,2 defined by Gw = VgE(w — W) is linear and 

also continuous by Corollary 4.4: 


||Gr 


"l2(c) 


< Clllw-Wll^: 


Wl/2,2 


<C2 


'^"^?.l/2.2 


and this implies that Gvk G4'(u) in T| 2 (c)j i-S-, VgE{vk — Vk) VgE(4'(u) — 'I'(u)) in T| 2 (c)- Now we 
show that indeed v = 'h(it). 


6.1 Identification of t; = 'I'(m) 

Let us define 


Jk{v) 


foImGkiv) HGkMeLl, 
0 otherwise 


and J(v) 


fofmGiv) 

0 


if Civ) € Ll, 
otherwise. 


Note that if n e L ^2 then Gk{v), G{v) G (see (43)). 

Lemma 6.2. The map 

V 1 -^ f f G{v) 

Jo Jr{t) 

from L ^2 into K is lower semicontinuous. 

Proof. First, observe that G: R —>■ K is convex, proper and continuous, hence (for example by adapting 
Proposition 8.1 in [49, Chapter II]) the map 

w !->■ / G{w) for w G L^(r(t)) 

Jr{t) 

(which is well-defined, for example, see (43)) is lower semicontinuous for each fixed t. li Vn ^ v in we 
have Vuj {t) —>■ v{t) in L^(F(t)) for almost all t, so 


f G(n(t)) < liminf [ G(n„^.(t)). 

Jr(t) rij-^oo 


(47) 


Integrating (47), and since Jp^^^ G{vn^{t)) > 0 and t !->■ Jp^^^ G(v„^.(<)) = /p^ G{vnj{t))J^ is measurable, we 
can apply Fatou’s lemma to give 

[ [ G{v{t)) < [ liminf [ G(n„j(t)) < liminf [ [ G{vn^{t)). 

Jo Jr{t) Jo Jr(t) rij^oo Jq Jp^jj 
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Thus far we have shown that for any sequence Vn ^ v converging in T(f) < liminfj_>oo holds for 

a subsequence rij. Now, if Vn ^ v in then it follows that there is a subsequence Vn^ such that 

liminf J(n„) = lira (48) 

n—^oo j—^oo 


by definition of the liminf (J is non-negative, so either liminf J(v„) = oo or liminf J(r;ji) = C > 0; the 
former case makes the problem trivial). We know that there is a subsequence of Uj such that J{v) < 
liminffe_>.oo J(vnj^) = limj_>.oo J(vnj) = liminf„_>.oo J(vn) with the first equality because the limit of 
is the same as the limit of J(vnj) and the second equality from (48). □ 

Lemma 6.3. We have u = 4'“^(r!). 

Proof. By convexity of Gk and G, Jk and J are also convex (see [10, §2.4]). If the Gateaux derivative of Jk 
or J exists at a particular point, then the set of subdifferentials of Jk or J coincides with the set of Gateaux 
derivatives at that point [42, Proposition 3.33]. By a direct calculation, the subdifferentials are 

dJk{vk) = {w G L \2 I w = 4'^^(vfc) in L]^^} and dJ{v) = {w G \ w = in LI 2 }. 


By definition (see [42, Definition 3.31]), since G dJk{vk), for all w G L| 2 , 


' y ' ' y ' ' ^ ' ' J I 

f [ Gk{vk)+ [ [ 'f’f^{vk)w< [ [ Gk{w)+ f ( '^f^(vk)vk. 

Jo Jr{t) Jo Jr(t) Jo Jr(t) Jo Jr(t) 


F 

/o JT{t) 


T 

0 Jr{t) 


F 

lo Jr(t) 


(49) 


We want to pass to the limit in this inequality using (46) and the methods of [30]. For the first term on the 
right hand side: for almost all t and almost all x G r(t), Gk{w{t,x)) —>■ G{w{t,x)) by the convergence of 
4 /[7^ —We also have by (43) \Gkiw{t,x))\ < C'(l'u;(<,a;)]^ -I- \w{t,x)\), and the right hand side is in L^i, 
so by the DOT, Gk{w) —>■ G{w) in L]^i, which obviously implies 


lo Jr(t) 


Gk{w) 



G{w). 


For the second term on the right hand side, since u G L‘j^ 2 , 



{'^k^M,Vk)L'^ 


7'2 


{u,v)i 


-1/2,2 ’-^^1/2,2 



For the first term on the left hand side, we first show an intermediary step, that 


lim [ [ Gk{vk) - G{vk) = 0. 

k^ooJo JT(t) 


(50) 


To see this, note that 


\Gk{vk{t,x)) - G{vk{t,x))\ 


Vk{t,x) 


(4/^(5)-4^-i(s)) 


< C sup |4'y(s) - 4- i(s)[, 
sG[-C.C] 


hence 


lo Jr{t) 


Gk{vk) - G{vk) 


< 1F|T sup l4'/\s)-4'-^(s)lC'^0. 

sG[-C.C] 


By weak lower semicontinuity of the map v 1 — > jJ" fr(t) (Lemma 6.2), we have 

/ / G(u) < liminf [ ( G(nfe) = liminf [ [ Gk{vk) 

Jo Jr(t) k^oo Jq Jr(t) k^oo Jq Jr(t) 
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with the equality by (50). Lastly, the second term on the left hand side is obvious. Now we can take lim inffc_>oo 
in (49) and use the above facts to get 






G{w) + 



uv, 


which is exactly the statement u S dJ{v), i.e., u = 4* ^(u). □ 

That u G follows from the strong convergence in ^ and the L°° estimate (42), see [2, §3(b)] for 

more details. Integrating (41) by parts over time and letting 77 G L^) with ri{T) = 0, the equation 

we want to pass to the limit in is 


/ / ii{t)uk{t) + / W g(^t)£t{'^k{uk{t)))W g(t)E{t)r]{t) = / uor]{0), 

/o Jr{t) Jo Jc(t) Jfo 


and this is easily done using the convergence results and will result in the equation in Definition 1.1. 


6.2 Contraction principle and conservation of mass 

We know that the solutions uik and U 2 k of the non-degenerate problem (with nonlinearity 4//;) and initial 
data uoi and uo 2 respectively satisfy 

/ {uik{t) - U 2 k{t))'^ < / (M 01 -W 02 )’'" for all fc (51) 

dr(t) Jto 

by Theorem 1.4. We have shown that (for a subsequence) utk converges to Ui, the solution of the fractional 
porous medium equation with initial data uoi. Now, with Uik := 4>-(^.)Uik, the bounds (42) and (45) translate 
into (see [2, Lemma 2.8] and [3, Proof of Theorem 2.33]) 

ll^jfc|lL°“(0,T;L“(ro)) + ll'“ifcllL2(o,T;lV-i/2,2(rj,)) < C, 

thus by Aubin-Lions (Theorem 11.5.16 in [15]), for a subsequence and for every t G [0,T], 

Uik{t) Ui{t) in W^/^’^(ro). (52) 

By the uniform bound, we have for almost all t that Uii^(k)it) — U 2 i-t(k)(t) ui(t) — U 2 {t) in L^(ro) (the 
identification is thanks to the strong convergence (52)). Since (•)^ is a convex function. It'. L^{r{t)) —?> K. 
defined by Itiw) = w'^ is convex, and clearly it is also continuous. Then, by a corollary of Mazur’s lemma 

[17, Corollary 3.8 and Remark 5], It is weakly lower semicontinuous, which from (51) gives 


/ (ui{t) — U 2 {t))'^ < / (uoi — uo 2 )^ for almost alH G [0,Tj. 

JT{t) Jto 

In fact this holds for every t G [0,T]. Take an arbitrary t G [0,T] and a sequence tg —>• t such that Uik{tj) 
is bounded in L°°(ro). This gives Uik{tj) ^ Uik{t) in L°°(ro) since Uik G (^^([O, T]; lT“^/^’^(ro)). The 
weak-star lower semicontinuity of norms gives 

||uife(i)|lLcx,(ro) < C for every t G [0,T], 

and the argument previously given can be repeated and we will get Uik{t) and then we can pass to 

the limit in the contraction result satisfied by uikitg) — U 2 k{tj), Hrst in j and then in k. 

The conservation of mass follows easily by passing to the limit in Jp^^^ Uk{t) = /p^ uq. 
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7 Concluding remarks 

The (non-fractional) porous medium equation on an evolving surface can be also tackled in this way, as a limit 
of approximations; of course the problem is easier in that case as we would not need §2, §4 and parts of §3, and 
the non-degenerate problem in §5 can be handled with a fixed point argument using the linear theory in [3], as 
done in [2] for a Stefan problem. We name a few of the many interesting open issues left. We required bounded 
initial data for the results above but the L^-continuous dependence result leaves us in good position to extend 
the results to integrable data if we manage to obtain a smoothing effect (for which the work [12] by Bonforte 
and Grillo may be useful). There is also the fast diffusion or the singular case where m G (0,1) which we 
have not addressed. A fundamental property enjoyed by solutions of the fractional porous medium equation 
on a stationary domain is regularity in time [26, Theorem 2.3], that is, the solution has a time derivative in 
L^. In the stationary case, this regularity is obtained partially by a rescaling argument of [8] and using the 
L^-continuous dependence applied to a solution and its rescaled version. This does not work in our setting 
since rescaled solutions live on a different evolving hypersurface, so the continuous dependence inequality 
cannot be applied. This result would be useful because it would allow us to study qualitative properties such 
as the effect the geometry of the hypersurface has on the solution. An obvious further extension is to study 
this theory of weak solutions with a general exponent in the fractional Laplacian (—App))'*: for this of course 
[19] is the obvious starting point and the methodology we used in this paper should work. 
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